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ABSTRACT 

We present a high-resolution, multi-wavelength study of the massive protostellar cluster NGC 6334 1(N) that 
combines new spectral line data from the Submillimeter Array (SMA) and VLA with a reanalysis of archival 
VLA continuum data, 2MASS and Spitzer images. As shown previously, the brightest 1.3 mm source SMAl 
contains substructure at subarcsecond resolution, and we report the first detection of SMAlb at 3.6 cm along 
with a new spatial component at 7 mm (SMAld). We find SMAl (aggregate of sources a, b, c, and d) and 
SMA4 to be comprised of free-free and dust components, while SMA6 shows only dust emission. Our 1 "5 
resolution 1.3 mm molecular line images reveal substantial hot-core line emission toward SMAl and to a lesser 
degree SMA2. We find CH3OH rotation temperatures of 165 ±9 K and 145 ± 12 K for SMAl and SMA2, 
respectively. We estimate a diameter of 1400 AU for the SMAl hot core emission, encompassing both SMAlb 
and SMAld, and speculate that these sources comprise a > 800 AU separation binary that may explain the 
previously-suggested precession of the outflow emanating from the SMAl region. Compact line emission 
from SMA4 is weak, and none is seen toward SMA6. The LSR velocities of SMAl, SMA2, and SMA4 ah 
differ by 1-2 km s"^ Outflow activity from SMAl, SMA2, SMA4, and SMA6 is observed in several molecules 
including SiO(5^) and IRAC 4.5 fim emission; 24 fim emission from SMA4 is also detected. Eleven water 
maser groups are detected, eight of which coincide with SMAl, SMA2, SMA4, and SMA6, while two others 
are associated with the Sandefl et al. (2000) source SM2. We also detect a total of 83 Class I CH3OH 44 GHz 
maser spots which likely result from the combined activity of many outflows. Our observations paint the 
portrait of multiple young hot cores in a protocluster prior to the stage where its members become visible in 
the near-infrared. 

Subject headings: stars: formation — infrared: stars — ISM: individual (NGC 6334 1(N)) — masers — tech- 
niques: interferometric 



1. INTRODUCTION 

J> ' The formation process of massive stars is an important but 
poorly-understood phenomenon. Because most regions that 
'/"J are forming massive stars lie at distances of several kilo- 
parsecs, and most massive stars form in clusters, the iden- 
tification of high mass protostars requires both good sensi- 
^\ tivity and high angular resolution to resolve the continuum 
and line emission within young protoclusters. The advent of 
^\ (sub)millimeter interferometers allows the study of deeply- 
embedded thermal gas emission at arcsecond angular reso- 
ld/ lution, a realm that has historically been accessible only via 
. !^ maser and cm A continuum emission. Research that com- 
1 bines observations of dust, thermal gas and maser gas offers a 
■ more complete picture of the star formation process, as these 
' phenomena potentially trace objects at different evolutionary 
stages. 

NGC 6334 is a luminous (> 10* Lq) and relatively nearby 
(1.7 kpc: ■NeckeL1978) region of massive sta r formation con- 
tainin g sites at various s t ages o f evolution (iStraw & Hvlandl 
119891) . [Matthews et al.1 (120081 ) have recently completed a 
comprehensive single dish study of the submillimeter dust 
properties of the whole complex with ^ 14" resolution and 
find a total dust mas s of ~ 1.7 x 10'* Mq. Recent work by 
iPersi & Tapial (120081) suggests fliat the distance to NGC 6334 
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may be as close as 1 .6 kpc, though we adopt 1 .7 kpc for ease 
of comparison with earlier work. The d ust core "I(N)" was 
first detect ed at 1 mm b y ICheung et al.l (11978.) and later at 
400 fj,m bv lGezaril (119821) toward the northern end of the com- 
plex. 

In the past decade, direct evidence of protostars in the I(N) 
region were found in the form of Class I and II methanol 
masers (Kogan & Slysh 1998; Walsh et al. 1998) and a faint 
3.6 cm source (Car ral et al.l |2002|) . In a previous paper, 
we presented the first millimeter interferometric observa- 
tions of I(N) with ^ 1 "5 resolution which revealed a clus- 
ter of compact dust continuum cores likely to be a proto- 
cluster (Hunter etal. 2006, see also Figure 1). Subsequent 
subarcsecond resolutio n 7 mm continuum observations by 
iRodrfguez et al.l (l2007h resolved the principal source SMAl 
into several components, and d etected possible co unterparts 
to SMA4 and SMA6 (also see ICarral et"an 12002!) . Our re- 
cent observations of the ammonia inversion transitions up to 
(6,6) demonstrated the presence of compact clumps of hot 
thermal gas within this protocluster (Beuther et al. 2007). Al- 
though water masers are excellent tracers of high mass star 
formation, previous obser vations of this region have b een lim- 
ited to single dish data (iMoran & Rodrigue3 1 198Cft . Like- 
wise, most of the millimeter line studies of this region have 
been perforrn e d with single-dish telescopes. For example, 
iPirogov et al.' ('2003') imaged the N2H+ (1-0) fine with a 
beamsize of 40" and found significantly stronger emission to- 
ward I(N) than toward the more luminous massive protoclus- 
ter (source I) that lies ^ 1' south of I(N). In addition, the SiO 
(2-1) and (5^) transitions show strong wing emission from 
-50 to H-40 km s"' when observed with beamsizes of 58" and 
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23", respectively (iMegeath & TieftrunJdl 19991) . 

Recently, we reported 3 mm interferometric imaging of 
I(N), which revealed compact CH3CN (5^) emission to- 
wa rd SMAl along with a bipolar outflow seen in HCN(1- 
0) (iBeuther et alj|2008l) . In this paper, we present the de- 
tails of our 1.3 mm Submillimeter Array^ (SMA) spectral 
line dataset, along with new NRAO Very Large Array*^ (VLA) 
H2O maser and 44 GHz Class I CH3OH maser data. We have 
also re-analyzed archival VLA continuum data from 3.6 cm to 
7 mm, 2MASS near-IR and Spitzer mid-IR IRAC and MIPS 
images. Together these data confirm the identification of 
NGC 6334 I(N) as an actively forming massive protostellar 
cluster. 

2. OBSERVATIONS 
2.1. Submillimeter Array (SMA) 

Our SMA observations were made with six antennas - in 
May 2004 in the compact configuration, and in May 2005 
in the extended configuration. In both tracks, two pointings 
were observed: source I at 17*'20'"53.44\ -35°47'02.2" and 
source I(N) at 17^20'"54.63^ -35°45'08.5" (J2000). Only 
the data for source I(N) are presented in this paper. The SMA 
receive rs are double sideba nd SIS mixers with 2 GHz band- 
width (iBlundell et alJlT998h . The local oscillator was tuned 
to provide 216.6-218.6 GHz in LSB and 226.6-228.6 GHz in 
USB, yielding a FWHP primary beam size of ^ 55". The cor- 
relator was configured for 3072 channels per sideband with a 
uniform spectral resolution of 0.8125 MHz. The correlator 
sideband separation was about 18 dB at the time of these ob- 
servations. Typical system temperatures were 200 K. The pro- 
jected baseline lengths ranged from 10 to 180 m. A continuum 
image produced from these data was presented in Hunt er et al.l 
(|2006), however in order to correct for a half-channel error 
in SMA velocity labeling discovered in November 2007, the 
data have been completely re-reduced. 

The gain calibrators were the quasars NRAO 530 (23" 
distant) and J1924-292 (27" distant). Bandpass calibration 
was performed with observations of Uranus (2004 data) and 
3C279 (2005 data). Flux calibration is based on observations 
of Jovian satellites and regular SMA monitoring of the quasar 
flux densities. The data were calibrated in Miriad, then ex- 
ported to AIPS for imaging; the two sidebands were reduced 
independently. The AIPS task UVLSF was used to sepa- 
rate the line and continuum emission. Self-calibration was 
derived from the continuum, and solutions were transferred 
to the spectral line data. After self-calibration, the contin- 
uum data from both sidebands were combined and imaged 
with both natural and uniform weighting in order to achieve 
the best surface brightness sensitivity (the former) and best 
angular resolution (the latter). The naturally weighted im- 
age has an angular resolution of 2"2 x 1"3 at PA. +15° and 
the uniformly weighted image has an angular resolution of 
1"8 X 0."9 at PA. +8°. The rms noise levels are 8.4 and 6.4 
mJy beam"' for the natural and uniform weighted images, re- 
spectively. The peak signal-to-noise of the naturally weighted 
continuum image has increased to 1 26 co mpared to 1 14 for 
the image presented in lHunter et al.l (120061) . The spectral line 

' The Submillimeter Airay (SMA) is a collaborative project between the 
Smithsonian Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy & Astrophysics of Taiwan. 

* The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under agreement by the Associated Universi- 
ties, Inc. 



data were resampled to a spectral resolution of 1 . 1 km s and 
imaged with natural weighting; the rms noise per channel is 
180 mJy beam"'. The absolute position uncertainty of these 
data is estimated to be 0"3. 

2.2. Very Large Array (VLA ) 
2.2.1. Water Maser Observations 

The 1.3 cm (22.235 GHz) H2O maser data were taken with 
the VLA on 2006 July 03 (project AH915) with onHne Han- 
ning smoothing and one polarization (RR). The full band- 
width of the data is 3.125 MHz (40 km s"' usable) and the 
channel separation and resolution is 24.4 kHz (0.33 km s"'). 
The absolute flux scale was set assuming a 22.235 GHz flux 
density of 2.54 Jy for 3C286. Fast-switching was employed 
using J 17 17-337 and a cycle time of 2 minutes. The calibrator 
J 1924-292 was used for bandpass calibration. The pointing 
model was updated once per hour during the observation. Af- 
ter flagging periods of poor phase coherence (due to weather), 
the on-source time is approximately 60 minutes. The FWHM 
primary beam of these data is 120" and the synthesized beam 
is0"79x0."25 atPA.=H-7°. 

It was necessary to remove the remarkably strong masers 
detected toward NGC 6334 I in the sidelobes (112" from the 
phase center of I(N)) in order to achieve high dynamic range 
on the target masers in source I(N); source I masers were de- 
tected with peak intensities (uncorrected for primary beam 
attenuation) as high as 112 Jy beam"' despite their location 
near the first null. The removal was accomplished by self- 
calibrating the data on the strongest source I channel, apply- 
ing this calibration to the full line dataset, imaging both fields, 
subtracting the clean components from source I (only) in the 
M-v plane (UVSUB), inverting the calibration (CLINV), re- 
self-calibrating on the strongest I(N) maser channel, applying 
these new corrections and making the final I(N) line cube (this 
process is often called "peeling"). This procedure improved 
the dynamic range in the I(N) channels at the same velocities 
as the strong source I masers by a factor of ~ 3. The rms noise 
in a single channel is 6 mJy beam"' . The absolute position un- 
certainty is approximately 0"01. Flux density measurements 
were made on an image corrected for primary beam attenua- 
tion. 

2.2.2. 44 GHz Methanol Maser Observations 

The 7 mm (44.06941 GHz) Class I CH3OH (7o,7-6i,6) A+ 
maser transition was observed on 2008 February 18 with the 
VLA in CnB -configuration (project code AC904). Correla- 
tor mode 2AB was used with a velocity channel width of 
0.17 km s"'. NGC6334I(N) was observed for 7.5 minutes 
(on-source) as a check source for a larger survey of massive 
protostellar candidates ('Cvga nowski et al.]|2009h . The gain, 
bandpass, and flux calibrators were J1717-337, J2253-I-161, 
and 3C286, respectively. Fast switching was employed with a 
cycle time of 2 minutes. The absolute positions are expected 
to be better than 0"2. The synthesized beam is l."30 x 0."40 
at PA. = -44°, and the velocity resolution is 0.33 km s"' af- 
ter Hanning smoothing. Due to the snapshot nature of these 
data, and consequent poor u-v coverage, the noise is dynamic 
range limited in channels with strong maser emission. Chan- 
nels with relatively weak maser emission have rms noise lev- 
els of ~ 100 mJy beam"', while it is up to 5 x poorer in chan- 
nels with strong maser emission (i.e. peak > 50 Jy). 

2.2.3. Archival 7 mm, L3 cm, and 3.6 cm Continuum Observations 
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Fig. 1. — Composite image showing the naturally weighted SMA 1.3 mm continuum emission from I(N) with surrounding insets showing the VLA H2O 
(22.235 GHz) maser spectra. The 1.3 mm continuum image has a resolution of 2"2 X lf'3 (P.A.=+15°) and the contour levels are 30 (3.5(t), 60, 120, 180, 240, 
320, 400, 640, 800, 960 mJy beam"' . The blue crosses show the intensity weighted centroid positions of the water maser groups listed in Table |2] (derived from 
Table[3). The water maser spectra show the average emission from each group. The dashed arc indicates the SMA primary beam, and the dashed arc indicates 
the VLA primary beam at 22 GHz. The 1.3 mm image has not been corrected for primary beam attenuation. 



We have re-reduced the archival VLA 7 mm and 
1.3 cm continu um observations originally presented by 
[Rodriguez et alj (i2007, project codes AZ159 and AZ152, fur- 
ther details of these observations can be found in that pa- 
per). We achieved a significantly better rms noise in the 
7 mm image (210 /iJy beam"') than those authors report (320 
^Jy beam"'). This improvement is likely due to the fact that 
we simultaneously imaged a field toward source I which con- 
tains a stron g ultracompact H n region (2' south of I(N), see 
for example [Hunter et al.ll2006 ): images made without clean- 
ing this Hn region show significant sidelobes at the loca- 
tion of I(N). We also used the MPS task VLANT to reti'ieve 
and apply post-observation antenna position corrections to the 
data, and used a model for the brightness distribution of the 
flux cahbrator 3C286. The flux densities derived for the phase 
calibrator (J 1720-355) are within 2% of those reported by 
[Rodriguez et al. (2007). The angular resolution of the 7 mm 
data is 0."65 x 0"44, P.A.=+25°. The 1.3 cm data were re- 
duced in a similar manner and the rms noise is 64 /iJy beam"' , 



similar to that obtained by [Rodriguez et all ( [2007h : the angu- 
lar resolution is 0"38 x 0."27, P.A.=+19°. The 1.3 cm and 
7 mm data show remarkable position agreement given that 
they were observed on different days - a testament to the ac- 
curacy of the VLA when reference pointing and fast switching 
are employed. We estimate that the absolute position accura- 
cies for these two datasets are better than 0"05. 

We have also re-reduced the archival VLA 3.6 cm c ontin- 
uum observations originally presented bv [Carral et alj ([20021 
project code AM495, further details of the observing parame- 
ters can be found in that paper). The pointing center for these 
data is located 1 .5' south of source I(N) (the 3.6 cm FWHP is 
5. '4). These data were re-reduced as described above for the 
7 mm and 1.3 cm data. In particular, a model for the bright- 
ness distribution of the flux calibrator (3C286) w as employed 
- a me thod not available at the time of the original [Carral et all 
([2002h reduction. The average flux densities (between the 
two IFs) derived for the two phase calibrators, J1733-130 and 
J1751-253 ai-e 6.67 ±0.03 and 0.267 ±0.001, respectively. 
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In contrast ICaral et alj ( l2002h obtained 7.6 ± 0.2 for J1733- 
130, 14% higher than our more accurate result (a flux density 
for J1751-253 was not reported by these authors). After self- 
calibration, the final image was made with baselines longer 
than lOOkA to minimize artifacts from the nearby (within 2') 
extended H II regions "I" and "E", which fall within the VLA 
3.6 cm FWHP primary beam. The resolution of the final 
3.6 cm image of the 1(N) region is 0"51 x 0."20, P.A.=-1° and 
the rms noise is 41 ^Jy beam"' after primary beam correction 
jCarral et al. 2002 report an rms noise of 60 ptJy beam"' but 
it is unclear whether they applied primary beam correction). 
The 3.6 cm data did not employ reference pointing or fast 
switching, and we estimate an absolute position accuracy for 
these data of 0"3. We note that there is a consistent 0."2 offset 
(to the SE) between the 3.6 cm data and the 1.3 cm and 7 mm 
data; we have not shifted the images to match. 

2.3. Archival Australia Telescope Compact Array (ATCA) 
3.4 mm Observations 

We have r e-imaged the 3. 4 mm ATCA continuum data 
reported by 'Beut her et al.l (120081) . The rms noise is 
3.3 mJy beam"' (after primary beam correction) and the res- 
olution is 2"3 X l."8, P.A.=83°. Due to the small primary 
beam of the ATCA at this wavelength (FWHP 33"), primary 
beam correction for this image is critical. This accounts for 
the apparent increase in rms noise compared to that reported 
by Beuther 'et al.l ( 120081 2.6 mJy beam"'; our rms noise is 
2.2 mJy beam"' before primary beam correction); the peak 
signal-to-noise has increased by a factor of 1.3. We estimate 
an absolute position accuracy for these data of 0."3. 

2.4. Archival Infrared Data 

For near-infrared p hotometry, we anal yzed archival 
2MASS' survey images ( iSkrutskie et al.ll2006 ^. We have also 
examined all of the Spitzer a rchival data for the 1(N) region. 
Of the available IRAC data (Fazio et al. 2004), the data ob- 
tained in project 46 on 2004, August 2 by PI Fazio are by far 
the deepest (666 s). In this paper we use the PB CD images 
produced by pipeline version 14.0 Hunter et al. (a different 
version of these data were also presented in 2006). Since sat- 
uration in the 1(N) region was not an issue, only the long inte- 
gration (10.4 s) data are used. Unfortunately, of the available 
MIPS data, only the 24 ^m data are usable - the 70 and 160 
fim data show tremendous striping and other artifacts. In this 
paper we utilize the 24 /x m data from the MIPSGAL Legacy 
survey (ICarev et al.ll2009h . 

We performed aperture photometry on all infrared images 
by defining irregular polygonal shapes (the "source" shape) 
surrounding each extended knot of emission (presumed to be 
shock fronts from one or more outflows, not point sources). A 
background is taken as the mode of pixel values in an scaled 
annulus, between 1.5 and 2 times radially expanded from the 
centroid of the source shape. These polygons are translated to 
each infrared image by their locations on the sky, so sample 
the identical regions even though the images have different 
pixel sizes and world coordinate systems. Fluxes were con- 
verted to Jy using the photometric information for each ob- 
servatory. 

' Two Micron All Sky Survey, a joint project of the University of Mas- 
sachusetts and the Infrared Processing and Analysis Center/California Insti- 
tute of Technology, funded by the National Aeronautics and Space Adminis- 
tration and the National Science Foundation 



2.5. Caltech Submillimeter Observatory ( CSO) 

Our CSO'*' observations were obtained on 13 April 1993 
using the facility 3 45 GHz SIS receiver with lead junctions 
dElhson etal.lll989h and a 500 MHz AOS backend. Several 
on-the-fly maps of the CO 3-2 line were obtained on a 10" 
grid, providing near-Nyquist sampling of the 20" beam of the 
telescope at this frequency. At the end of each row, a desig- 
nated off position (17''24™56.06^ -36°21'40.1" J2000) rela- 
tively free of CO emission was observed to provide the spec- 
tral baseline. After combining the maps, the total area cov- 
ered was 4' X 4' (25 x 25 grid), with 30 seconds of on-source 
integration per point in the central portion containing source 
1 and 1(N), and 15 seconds per point on the periphery. The 
DSB receiver temperature was measured to be 320 K using 
the y-factor method, and the typical SSB system temperature 
was 1800 K at an elevation of 30°. The spectra have been cor- 
rected to the main beam brightness temperature scale using an 
efficiency of 68% as measured on Jupiter on the same night. 

3. RESULTS 

Our multiwavelength images reveal several sites of star for- 
mation activity within the 1(N) region. In the following sec- 
tions, we describe the variety of observed phenomena. In the 
following sections linear offsets are given with the > sym- 
bol to indicate that the measured values are lower limits since 
only the projected separation is available (a distance of 1.7 
kpc is assumed). 

3.1. SMA 1.3 mm Continuum 

The naturally weighted 1.3 mm SMA continuum image of 
I(N) is presented in Figure [T] (cent er). This image is quali- 
tatively similar to that presented in Hu nter et al.l (12006 ). but 
the calibration has been refined for improved signal-to-noise. 
These data resolve the presence of at least seven 1. 3 mm 
"cores" in I(N). As pointed out in iHunter et al.l (12006; ). this 
region has the hallmarks of a "protocluster", i.e. a region that 
will form a cluster of intermediate to massive stars. The in- 
tegrated flux density within the 3<t contour (30 mJy beam"', 
see Fig. [Til of the naturally weighted 1.3 mm continuum im- 
age is 7.4 ± 0.8 Jy (this is 1 .6 times higher than reported in 
IHunter et al.ll2006h ; this estima te includ es both statistical and 
10% calibration uncertainties. ISandelll (12000) find a 1.3 mm 
single dish flux density of 14.2 ± 0.2 Jy, implying that we have 
recovered about 50% of the total flux with the SMA. 

We have used the uniformly weighted image corrected for 
primary beam attenuation to fit the seven peaks identified in 
Figure[T]with a single 2-D Gaussian (using JMFIT in AIPS). 
The fitted positions and flux densities of the 1.3 mm cores 
are Hsted in Tabled SMAl, SMA4, and SMA6 in particular, 
are not well-fit by a single Gaussian component, and instead 
appear to have a compact (unresolved, FWHM < 1") region 
of emission, embedded in a more extended though still com- 
pact (FWHM < 3") region. However, attempts to uniquely 
disentangle these two (or more) components with the current 
resolution and sensitivity were not successful. The JMFIT 
sizes listed in Table [T] are representative of the more extended 
(though still compact on the scale of the I(N) protocluster) 
component, and the integrated intensities and other quantities 
calculated from it (7/, and mass) are likely to be underesti- 
mates. These sizes while quite significant from a signal-to- 
noise perspective should be viewed as upper limits. 

The Caltech Submillimeter Observatory 10.4 m is operated by Caltech 
under a contract from the National Science Foundation. 
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Fig. 2. — Zoomed multiwavelength images of I(N). (a) Zoom shows the SMAl, SMA2, and SMA3 regions with thick contours showing the uniformly weighted 
SMA 1.3 mm continuum image with contour levels of 170, 250, 330, 490, and 650 mJy beam"' . Thin contours and greyscale show the VLA 7 mm continuum 
emission (contour levels are -0.6, 0.6, 0.8, 1.0, 1.2 mJy beam"'). The water maser positions are shown as colored crosses with a con'esponding color coded 
velocity legend; note that the color scale was set using an assumed -3 km s"' systemic velocity for SMAl. (b) Zoom of the SMAl region at 3.6 cm (dashed 
contours), 1 .3 cm (greyscale), and 7 mm (solid contours). The 3.6 cm contour levels are 0. 12, 0.20, 0.28, and 0.36 mJy beam"' , and the 7 mm contour levels are 
same as in (a). The water maser positions are also indicated with the same velocity color coding as (a). The position of the 44 GHz CH3OH maser is indicated 
by the yellow ellipse and it has a velocity of —0.9 km s"' . For both types of maser symbol size indicates position uncertainty. 



3.2. ATCA 3.4 mm Continuum 

As first discussed bv lBeuther et al.l ( 120081) . the 3.4 mm con- 
tinuum emission toward I(N) (not shown) is qualitatively sim- 
ilar to that at 1.3 mm . With our new imaging of the original 
iBeuther et al.l ([2008") data, we have been able to increase the 
signal-to-noise and reliably detect SMA4 at this wavelength. 
The fitted (JMFIT) parameters of the compact 3 mm contin- 
uum emission are given in Table [T] 

3.3. VLA Water Masers and Continuum Emission 

The association of water masers with molecular out- 
flows make them a good tr acer of protostellar activity 
(iTofani et al.l [T99l iHofner & Churchwell 1996) . We find 
strong 22 GHz H2O maser emission toward 11 distinct re- 
gions in NGC 6334 I(N), listed in Table|2]in order of increas- 
ing declination. Masers that lie in close proximity to each 
other (< 1") are collectively called by "group" names CI, 
C2...C11. The intensity weighted positions of these groups 
are listed in Table |2] and are plotted on the 1.3 mm contin- 
uum image in Figure [T] along with average spectra for each 
group (the position of each distinct H2O maser and its ve- 
locity are presented in Table [3] available in its entirety on- 
line). Eight of the groups lie within 2" of the 1.3 mm sources 
SMAl, SMA2, SMA4, and SMA6, whose positions are Usted 
in Table [T] Two of the maser groups (CI and C2) lie near 
the SW edge of the SMA primary beam. These masers are 
coincident with the single-dish millimeter continuum source 
SM2, identifie d in JCMT submillimeter continuum images by 
iSandelll (12000'). Maser group C4 appears in a region with no 
compact millimeter continuum emission at the present sensi- 
tivity. Instead it lies at the southern edge o f the 4.5 iim nebula 
that extends to the southwest from SMA4 dHunter et al.ll2006i 
also see §3.6). More detailed phenomenology for the SMA 
1.3 mm cores (except SMA5 and SMA7 which are weak and 
fairly diffuse) are given below. 



3.3.1. SMAl, SMAl, and SMA3 Regions 

Figure 12^ shows a close-up view of the region from SMAl 
through SMA3 in both 1.3 and 7 mm continuum. This region 
also contains H2O maser groups C6 through Cll. Follow- 
ing Rodriguez et al. (2007) we call lower frequency contin- 
uum emission that appears to be associated with a particular 
SMA source by the "SMA name". If more than one source 
is associated with a 1.3 mm source, the SMA name is ap- 
pended by a,b,c etc. Lower frequency continuum emission 
that is not clearly associated with any particular SMA source 
is called "VLA" followed by the number of the closest SMA 
source. From the improvement in t he 7 m m image quality 
over that shown in iRodrfguez et al.l (120071) . several new fea- 
tures are present: (1) The 7 mm source SMAlb is resolved 
into two distinct sources separated by ~ 0."45 (> 800 AU, 
also see Figure|2^), we shall refer to the eastern component as 
SMA Id; (2) SMA la is elongated NE/SW and appears to point 
toward SMAlb; (3) Two new 7 mm features are detected to 
the north of SMA2 and SMA3, the stronger of the two is clos- 
est to SMA3 (-- 1" north), and is thus called VLA3, while 
the other is barely above the 3cr level and is not discussed fur- 
ther The flux densities and positions of all of the continuum 
sources from 1.3 mm to 3.6 cm are provided in Table[T] 

Figure ^ further zooms in on the continuum emission 
from SMAlb and SMAld and the maser groups C8, C9, and 
CIO. From this figure it is clear that SMAld is also distinct 
from SMAlb at 1.3 cm. From our reanalysis of the 3.6 cm 
data originally reported by lCarral et al.l (12002) we also detect 
SMAlb at this frequency. As mentioned in §2.2 there appears 
to be a systematic 0"2 position offset of the 3.6 cm data which 
if removed would cause it to match the 1.3 cm and 7 mm po- 
sition for SMAlb. Maser group C8, located ~ 0."3 south of 
SMAlb, contains the brightest H2O maser in the whole region 
and is resolved into a linear north/south structure > 800 AU in 
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Fig. 3. — (a) Zoom of the SMA4 region with the grayscale and dashed contours showing the 1.3 cm continuum emission, dotted contours showing 3.6 cm, and 
the solid contours showing the uniformly weighted 1.3 mm continuum. The 3.6 cm contours are 0.12, 0.16, and 0.20 mJy beam"', the 1.3 cm contour levels are 
0.19, 0.25, 0.31, and 0.37 mJy beam"' , and the 1.3 mm contour levels ai'e 100, 120, 140, and 160 mJy beam"' . (b) Zoom of the SMA6 region with the greyscale 
and dashed contours showing 7 mm continuum emission and the solid contours showing the 1.3 mm continuum. The 7 mm contour levels are 0.6, 0.8, 1.0, 
and 1.2 mJy beam"', and the uniformly weighted 1.3 mm contour levels are 100. 140, 180, 220, and 260 mJy beam"'. The cross symbols show the C5 (a) and 
C3 (b) wate r maser positions with S/N ratio > 7, color-coded by velocity. The filled square symbols in (a) show the location of 6.7 GHz CH3OH masers from 
IWalsh et"an (r998), with same velocity coding as H2O masers. 



length. The velocity structure of this linear group of masers 
is broad and complex, although the northern part is mostly 
blueshifted while the southern part is mostly redshifted with 
respect to the ~ -3 km s"' systemic velocity of SMAl (see 
§3.5 for discussion of systemic velocities). These character- 
istics are similar to the brightest H7O maser c omponent (C6) 
in the protocluster AFGL 5142 dHunter et al. 'lggg) but on a 
somewhat larger physical scale. A 44 GHz maser with a ve- 
locity of ~ -0.9 km s"' is also coincident with this linear H2O 
maser feature. 

The redshifted and kinetically narrow H2O maser group C9 
is coincident with SMAld and contains the second brightest 
maser feature in I(N). A weak blueshifted group of masers 
(CIO) lies > 1200 AU southeast of SMAlb. Back on the 
larger scale of Fig. |2^, group CI 1 is another blueshifted nar- 
row feature that lies > 1300 AU northwest of SMAlc. Group 
C6 is > 400 AU northwest of SMA2 and is blueshifted with 
respect to the systemic velocity of SMA2 (^ -5 km s"'; see 
§3.5). Group C7 is a narrow but fairly strong component 
> 1500 AU northwest of SMA2, and is near the systemic ve- 
locity. 

3.3.2. SMA4 Region 

Looking further south toward SMA4, the spatial/velocity 
structure of water maser group C5, as well as the 1.3 mm and 
1.3 and 3.6 cm continuum emission in this region is shown 
in Figure [3^. There are three distinct spatial clusters of H2O 
masers within the C5 group, corresponding to the three dif- 
ferent velocity components. These three clusters of masers 
are oriented roughly north/south (P.A.=17°) with blueshifted 
emission to the north, and redshifted emission to the south. 
The C5 group Hes within 0."4 (> 700 AU) of the 1.3 mm to 
3.6 cm continuum emission peaks, just within the combined 



astronomical uncertainties of the data. The C5 group is also 
coincident with a cluster of Class 11 6.7 GHz CH3OH masers 
(Walsh et al. 1998; Caswell 200 9) which have a similar veloc- 
ity gradient. The 1.3 and 3.6 cm em issio n has been reporte d 
previously by lRodriguez et alJ (120071) and lCarral et"an (l2002h . 
respectively. We report the detection of SMA4 at 7 mm (not 
shown, see Table[T]i- As mentioned in §2.2 there appears to be 
a systematic offset of 0."2 of the 3.6 cm data which if removed 
would cause it to match the 1.3 cm and 7 mm position. It is 
currently unclear if the maser and 1.3 mm offsets are signifi- 
cant (or simply due to absolute position uncertainties). 

3.3.3. SMA6 Region 

A zoomed in view of the SMA6 region is shown in Fig- 
ure [3J). The C3 group of water masers is coincident with 
SMA6 and also has three distinct velocity clusters (see 
Fig. [3}'). The kinematics of C3 are curious, with the three 
clusters lying along a roughly north-south line with systemic 
velocity to the south, then redshifted in the middle, and then 
blueshifted emission to the north. The 7 mm continuum peak 
toward SMA6 lies > 700 AU north of the 1 .3 mm peak, co- 
incident with the northern blueshifted cluster of C3 masers. 
SMA6 is not detected at wavelengths longer than 7 mm. The 
1.3 mm morphology of SMA6 is suggestive of unresolved 
substructure in the north-south direction, and as mentioned 
in §3.1 is not well fit by a single Gaussian. Interestingly, al- 
though the 7 mm data has significantly better spatial resolu- 
tion than the 1.3 mm data, both images show a similar mor- 
phology. 

3.4. VLA 44 GHz Methanol masers 

Strong 44 GHz CH3OH maser emission from 
NGC 6334 I(N) was first detected in the single dish 
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Fig. 4. — SMA 1.3 mm continuum contours of I(N) (thin black contours) 
with thick black contours showing SCUBA 450 /im emission (G. Sandell, 
private communication). The 1.3 mm continuum image and contour levels 
are the same as Fig.[T] The SCUBA 450 fim resolution is 8" and the contour 
levels are 60, 80, 10 0, 120 16 Jy beam"' . The tip of the arrow for SM2 is the 
position reported by Sandell f3000). The colored crosses show the locations 
of 44 GHz Class I methanol masers from this work. The locations of the 
water maser groups (also see Fig.[T) are shown by I symbols and the 25 GHz 
CH3OH masers reported by Beuther et al. (2005) are shown by ♦ symbols. 
The X-ray sources from Feigelson et al (^009) are marked by magenta X's. 
The dashed arc indicates the SMA primary beam and the dotted arc indicates 
the VLA 44 GHz primary beam. The 1.3 mm SMA continuum resolution is 
shown in the lower left. 

survey of iHaschick et all (46" beam; I1990I) . iKogan & Slvshl 
(119981) resolved this emission into 23 spots using ten antennas 
of the VLA, but they were unable to determine accurate 
absolute positions. Our new observations with the full 
VLA provide a deeper more accurate view of the maser 
emission. By examining and comparing the details of the 
emission in each channel and from channel to channel, we 
have identified 83 individual maser spots. These are listed 
in Table ID available in its entirety online. The values of 
Vpeak refer to the center velocity of the channel in which the 
peak emission occurs. The tabulated positions are the fitted 
position in the peak channel for each component. The peak 
brightness temperatures of these masers (see Table |4|i lie far 
above the energy level of this transition (65 K, 43.7 cm"'), 
with values as high as 7], = 3.6 x 10^ K with the current 
spatial resolution. While we have detected far more spots 
of emissio n, the spatial and velocity distribution of the 
iKo^an &~Sl vsh (199 8|) spots are quite co nsistent with our 
data if one shifts the Kogan & Slvshl (Il998h positions by 7.3" 
to the south-southwest. At the current level of sensitivity, no 
thermal emission was detected. 

The maser positions are shown as velocity color-coded 
crosses in Figure 14.31 The peak maser velocities are in the 
range -8.49 to -0.35 km s~', close to the systemic velocity 
of the region. Figure |43] also shows a 450 /xm JCMT con- 
tinuum image with 8" resolution from I Sandelll (120001) . The 
masers are distributed over much of the area of dust emission, 
as delineated by the lowest 450 /im contour The masers can 
be divided into two main concentrations: the general area of 
the compact SMA continuum sources, a nd the area around 
the single-dish source SM2 (.Sandellll200 0). The former con- 
centration contains a more diverse range of velocities, though 
no obvious trends in kinematics versus position are seen. It 
is notable that while these masers are clustered around I(N) 
L3 mm continuum sources, only one is positionally coinci- 
dent with a compact SMA source (SMAl, see Fig. |2]i. The 
latter concentration towards the single dish continuum source 
SM2 is not associated with any 1.3 mm SMA continuum 
emission but this region lies outside the half-power point of 



the SMA data, significantly reducing the sensitivity in this re- 
gion. Two 25 G Hz Class I CH 3 OH m asers were also detected 
toward I(N) by i Beuther etall ( 12005 ). and both of these are 
coincident with 44 GHz maser emission and have comparable 
velocities as shown in Fig. 14.31 

The positi on of the X-ray so urces recently reported for this 
region (Fei gelson et al J 120091) are also shown in Figure |43] 
The median photon energy of many of these X-ray sources 
is relatively hard (E > 4 keV), suggesting deeply embedded 
sources. These authors have already noted the possible as- 
sociation of X-ray source 1454 with SMA7, and 1447 with 
SMA6 (although the latter is less likely given that its me- 
dian photon energy is only 1 keV). With our new data, we 
note four additional possible associations: X-ray source 1470 
is within 1"8 of the water maser component C2; and X-ray 
sources 1463, 1477 and 1487 each reside within a few arcsec- 
onds of a cluster of methanol masers. Source 1463 lies SSE 
of SMAl, just beyond the lowest millimeter contour Source 
1477 aligns with the shoulder of 450 /im emission northeast 
of SM2, while source 1487 lies just east of the lowest 450 yum 
contour 

3.5. SMA 1.3 Millimeter Molecular Line Data 

A total of 79 lines from 19 molecular species are detected 
toward SMAl (Table |5]l along with 14 unidentified features 
within the 4 GHz of total bandwidth observed by the SMA. 
The spectra toward the SMAl 1.3 mm continuum peak are 
shown in Figur e IH The lines were indentified using the 
CDMS (MuUer etd] l200l . JPL (Pickett et al. 1998), and 
Splatalogue" spectral line databases. The large number of 
organic transitions detected, coupled with modest line bright- 
ness temperatures of 3 K suggests that this source harbors 
a moderate "hot core" (Ivan Dishoeck & Blakelll998l) . Many 
of the stronger lines detected toward SMAl are also detected 
toward SMA2, though most are weaker, suggesting that this 
source also harbors a weak hot core. Other sources ob- 
served with similar linear resolution with this SMA setup have 
yielded both less diverse (S255N, Cy ganowski et al. 2007), 
as well as much richer hot core spectra (NGC6334I, Brogan 
et al., in prep). The broadband SEST single-dish molecular 
line survey toward NGC 63341 and I(N) also shows that while 
source I(N) exhibits significant hot core emission, source I 
shows a much richer a rray of molecular line species along 
with stronger emission dThorwirth et al.ll2003h . 

3.5.1. Extended Molecular Line Emission 

Many of the lines from low excitation states (£/ < 70 cm"') 
show emission extending to the NE and SW of the SMA con- 
tinuum emission (see Figure|6]l. This emission is clearly trac- 
ing one or more outflows emanating from the protocluster. 
Interpretation of the extended emission is complicated by the 
fact that the SMA data are not sensitive to smooth emission 
on sizescales > 20" (due to spatial filtering by the interfer- 
omter). Generally, the emission to the northeast is red-shifted, 
while that to the southwest is blueshifted, but the kinematics 
become very confused in the southern part of I(N) with both 
red and blueshifted emission present. SiO in particular shows 
a striking bipolar appearance centered on SMAl. Channel 
maps of SiO (Figure |7^) reveal the complex velocity field of 
the extended emission. A lobe of molecular material also ex- 
tends westward from SMA4 at low, but blueshifted veloci- 

' ' http://www.cv.nrao.edu/php/aremijan/splat_beta/ 
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ties in SiO and several other species shown in Fig.|6](H2CO, 
CH3OH, and H2S for example). 

A more compact SiO velocity gradient also appears to be 
centered on SMA6. CSO CO(3-2) spectra with 20" reso- 
lution show that the broad line wings toward SMA6 are of 
similar magnitude to those seen toward SMAl (see Fig. |7j)), 
suggesting both objects harbor outflows. The integrated red 
and blue-shifted emission around SMA6 is shown in detail in 
Figure |7};, and is more or less centered on the CN and H2CO 



(El = 7 cm"') absorption velocity of ~ -4 km s"' (see §3.5.2). 
This figure also shows that the sense of the SiO gradient is 
reflected in the position and velocity of the 44 GHz methanol 
masers surrounding SMA6. However, on smaller scales the 
C3 group of water masers is oriented in a more north-south 
direction (Fig. It is unclear how the curious kinematics 
of the C3 water masers (from south to north: systemic, then 
red, then blueshifted features) fits in with the compact NE-S W 
outflow. 
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Fig. 6. — In each panel, the colorscale is the 1.3 mm continuum and the 
contours show the integrated intensity of the indicated spectral lines. The 
dashed circle is the FWHP of the SMA primary beam; these images have 
not been corrected for primary beam attenuation. The contour levels are 
H2CO(3o.3-2o,2): 4.35,7.25, 10.15, 13.05 Jy km s"'; DCN: 1.77, 2.95,4.13 
Jy km s-l ; SiO: 6, 10, 14, 18, 22 Jy km s"'; CH3OH-E (4+2,2-3+1.2): 4.8, 
8, 11.2, 14.4, 17.6Jy kms-'; H2CO (32,2-22.i): 3.24, 5.4, 7.56 Jy km s"' : 
H2S: 2.55, 4.25, 5.95, 9.35 Jy km s"'. The number below each molecular 
name is the Slower of the transition in units of cm"' . 

3.5.2. Compact Molecular Line Emission and Absorption 

As seen in Figure [8] emission from the higher excitation 
lines (El > 70 cm"') is concentrated toward SMAl. In sev- 
eral species, a secondary peak of emission appears toward 
SMA2. Indeed, the SMA2 peak in methyl formate and the 
lower lying transitions of methanol is almost equal to that of 
SMAl. A few vibrational states of CH3OH and HC3N ai-e de- 
tected exclusively toward SMAl, suggesting that it is warmer 
and/or denser than SMA2. Several species also show weak 
emission toward SMA4 and to the north and west of SMA6. 
Figure |9] shows the first moment maps for several representa- 
tive species detected in emission, along with the CN molecule 
which is detected in absorption. Because of the possibility 
of confusion with outflowing or infalling gas, it is difficult to 
conclusively determine the systemic velocities of all of the 
1.3 mm cores. Both SMAl and SMA2 show enough con- 
sistency in a range of species with compact emission to esti- 
mate that y/,,(SMAl) = -2.8 ±0.5 km s"' and y/,,(SMA2) = 
-4.5 ±0.5 km s"'. Compact emission from SMA4 is weak 
and shows a peak emission velocity of V/sr(SMA4) = -5 ± 1 
km s~'. No emission centered on SMA6 was detected. 

The CN(2-1) lines are seen in absorption against most of 
the stronger compact continuum sources, while most of its ex- 
tended emission is resolved out by the interferometer None 
of the 1.3 mm continuum brightness temperatures (Th) are 
high enough at the current angular resolution to detect the 
El = 40 cm"' H2CO transition in absorption (negative) against 
it, though self-absorption is seen toward SMAl and SMA2. 
Only SMAl has a continuum Tt, (~ 15 K at current resolu- 



tion) high enough to detect the H2CO Ei = 1 cm ' transition 
in (negative) absorption. Figure [TOh and b show spectra of 
HCOOCH3 {El = 77 cm"'), H2CO at = 7 and 40 cm"', and 
CN(£'/ = 3.8 cm"') toward the SMAl andSMA2 1.3 mm con- 
tinuum peaks. The absorption velocity of CN toward SMAl 
is in good agreement with V/jr(SMAl) = -2.8 ±0.5 km s~'. 
The absorption seen in the H2CO (£/ = 7 cm"') line is slightly 
redshifted by 0.5 km s"' compared to CN and the V/j,.. Inter- 
estingly, the El = 40 cm"' H2CO line shows a self-absorption 
dip redshifted by ~ 1 km s"' relative to CN and the Visr (see 

SMA2 shows CN absorption redshifted by 4 km s ' from 
the V/j,.(SMA2) = -4.5 ±0.5 km s"', suggesting strong in- 
fall is present towards this source (see Fig. [TOb). In contrast, 
the El = 40 cm"' H2CO emission peak is slightly blueshifted 
('^ 0.5 km s"') compared to other emission lines (see for ex- 
ample HCOOCH3 in FiglH?), while the £/ = 7 cm"' H2CO 
transition's peak is blueshifted by 2 km s"'. Neither H2CO 
line is Gaussian in shape, with the redshifted side of both 
transitions cutting off very sharply as might be expected from 
very strong, red-shifted self-absorption. Weak blueshifted (-5 
to -15 km s"') H2O masers were also detected toward SMA2 
so these masers and H2CO could be tracing a pole-on outflow, 
with the Visr emission missing due to self-absorption and the 
redshifted H2CO outflow emission obscured by the contin- 
uum. 

No compact molecular line emission is detected toward 
SMA3 so it is unclear if the CN absorption toward this core 
with a velocity of ^ -2 km s"' is also tracing infall or the Visr- 
Toward SMA4, CN shows absorption at -5 km s"' in good 
agreement with the weak emission lines detected toward this 
source including £, = 40 cm"' H2CO (y;„(SMA4) = -5 ± 1 
km s"'). The £/ = 7 cm"' H2CO transition toward SMA4 
is difficult to interpret as it is non-Gaussian in shape with a 
peak at ^ -2.5 km s"' and a significant blue-shifted shoul- 
der The CN absorption toward SMA6 has a velocity of 
-4.0 ±0.5 km s"', in good agreement with the velocity of a 
self-absorption dip seen in £/ = 7 cm"' H2CO. Since no com- 
pact emission is detected towards this source it is unclear how 
this absorption velocity compares to the Visr of SMA6. 

Also of interest is the barely resolved SE-NW (or some- 
times more E-W) velocity gradient across SMAl evident in a 
number of compact species (see for example HCOOCH3 and 
CH3OH in Fig. |9]l. The full width of the velocity gradient 
is about ^ 5 km s"' (from inspecting the line cubes), i n rea- 
sonabl e agreement with the gradient inferred by Beuth er et alJ 
(l2007h for NH3(6,6) based on a double peaked profile. This 
gradient is discussed further in §4.5. 

3.6. Mid-infrared Emission 

Figure nn shows a three color mid-IR image of the I(N) re- 
gion with methanol maser positions and integrated SiO and 
1.3 mm continuum contours superposed. The extended 4.5 
lim emission to the SW of the I( N) continuum has been re- 
ported previously bv lHunter et al.l (|2006). In addition, a weak 
region of predominantly 4.5 fim emission is also visible be- 
tween the SMA6 continuum peak and the arc of 44 GHz 
CH3OH masers 5" to the SW; this region is coincident with 
the blueshifted side of the SMA6 outflow (Fig. |7j;). Clumpy 
knots of extended 4.5 /im emission are also present to the NE 
of SMAl, and in the vicinity of SMA7. We also present for 
the first time, an unresolved source of 24 /zm emission located 
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Fig. 7. — (a) Channel maps of SiO 5^. In each panel, the greyscale is the 1.3 mm continuum emission and the contours are SiO emission (levels = 0.29, 0.48, 
0.67, 0.86, 1.05 Jy). The outflow axis is drawn at position angle +39°. The dashed circle is the SMA primary beam, (b) CSO CO (3-2) spectra taken near SMAl, 
SMA6 and a position offset —40", +53" to the northwest (away from any obvious outflow emission). The SMAl spectrum was taken 2.9" from SMAl, and the 
SMA6 spectrum is the average of the three spectra taken closest to SMA6 (5", 7", and 8" away), (c) Close up of SiO (5^) outflow around SMA6. The contours 
show SiO (5^) integrated emission, red contours show integrated intensity from +17 to -6 km s"' and blue contours show integrated intensity from from -8 to 
-31 km s"' . The CH3OH and H2O maser positions are also shown. The solid black ellipse is the SMA synthesized beam. 

2" SW of SMA4. This figure demonstrates that with the 
exception of extended 4.5 /im emission, and a single source 
of 24 fim emission, NGC 6334 I(N) is relatively dim in the 
mid-lR. Indeed, SMA4 appears to be the only millimeter con- 
tinuum source directly associated with a mid-lR source. In 
contrast, the bright haze of red emission in the SW corner 
of Fig. [TT]emanates from the periphery of the saturated mid- 
lR bright protocluster NGC 6334 1. The nature of the mid-lR 
emission is discussed in detail in §4. 1 and §4.2. 

4. DISCUSSION 

4.1. Presence of Multiple Outflows 

Due to the ubiquity of bipolar outflows from protostars, 
one would expect to find multiple outflows in a protostellar 
cluster As discussed in Cyganows ki et al.l (12008) (and refer- 
ences therein), the 4.5 /im IRAC passband can be dominated 
by molecular line emission from H2 and vibrationally-excited 
CO in regions of strong shocks such as those found in mas- 
sive molecular outflows. To the NE of SMAl, the extended 
4.5 fim emission appears to be an extension of the large-scale 
NE-SW SMAl outflow ti'aced by S iO and other tracers of the 
extended emission (see Figs. |6|7bJTn i. To the south where 
the outflows from SMAl and SMA4 (westward flow) over- 
lap, it is unclear which outflow the 4.5 /im emission is aris- 
ing from. Interestingly, the 4.5 fim emission extends further 
than the SiO does to the west and SW. As described in §3.5, 
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Fig. 8. — In each panel, the greyscale is the 1.3 mm continuum and the 
contours are drawn from moment-zero images of spectral lines. The contour 
levels are: HCOOCH3: 2.1, 4.2, 6.3, 8.4, 10.5 Jy km s"'; HC3N: 2.55, 4.25, 
5.95, 9.35, 14.45, 19.55 Jy km s"'; C2H5CN: 1.8, 3.0, 4.2, 5.4 Jy km s"'; 
CH3OCH3: 1.35, 2.25, 3.15 Jy km s"'; CH3OH(220): 1.89, 3.78, 5.67, 
7.56 Jy km s"'; 2.31, 3.85, 5.39, 6.93 Jy km s"'; HC3N(v7=l): 1.5, 2.5, 
3.5 Jy km s"'; CH3OH(380): 2.1, 3.5, 4.9 Jy km s"'. The numbers below 
each molecular name is the iJiowcr of the transition in units of cm"' . 
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Right Ascension (J2000) 

Fig. 9. — In each panel, the color scale shows the first moment map of the 
indicated molecular line, while the solid black contours depict the integrated 
intensity image of the line emission. In the case of CN, the dashed contours 
indicate that the line appears in absorption. The levels are: HCOOCH3= 2.1, 
6.3, 10.5 Jy km s"'; HC3N= 2.55, 5.95, 14.45 Jy km s"'; CH30H= 1.89, 
5.67 Jy km s"' ; CN = -1.2, -2.4, -3.6 Jy km s"' . In each panel, the thin dotted 
contours show the 1.3 mm continuum emission (levels = 60, 180, 390, and 
900 mjy beam"'). The numbers below each molecular name is the Slower of 
the transition in units of cm"' . 
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Fig. 10. — (a) Sample spectra toward the 1.3 mm continuum peak of SMAl. 
(b) Sample spectra toward the 1.3 mm continuum peak of SMA2. In both 
plots the dashed line indicates the Vi,„ and the numbers in parenthesis under 
the molecule names are the lower state energy in cm"' . In (a) the two dotted 
lines indicate the velocity of the peak red and blueshifted H2CO (£; = 7 cm"' ) 
emission. In (b) the dotted line indicates the velocity of the peak blueshifted 
H2CO (El = 7 cm"') emission and the dot-dashed line indicates the velocity 
of maximum CN absorption. 

SMA2 may harbor a more or less pole-on outflow as traced 
by blueshifted H2O masers and possibly H2CO. As described 
in §3.6 and §3.7, a smaller scale NE-SW outflow also appears 
to emanate from SMA6, with blueshifted SiO and 44 GHz 
CH3OH maser emission located to the south coincident with 
4.5 /xm emission (Figs.lTbJTTTl. 



The ATCA molecular line data reported by iBeuther et alj 
( l2007i I2OO8I) (HCN(l-O) and NH3 for example) show good 
overall morphological agreement with the largescale emis- 
sion emanating fr om SM Al observed with the SMA (see 
for example Figs. I6I7I1 lb . In particular, the complex kine- 
matics of the extended emission were also noted by Beuther 
et al. and ascribed to possible precession. This idea is in 
good agreement with the discovery of a possible tight bi- 
nary in SMAlb and SMAld (separated by > 800 AU) if 
these two sources do in fact represent two different proto- 
stars (see^for example Chandler et al. 2005; Angladaetal] 
lIOOTt ICunningham et al.l 12009). Although Figs. il7}i,b,cIII] 
show SiO an d othe r mole cules in the vicinity of SMA6, 
IBeuther et all (120071 l2008h did not detect any significant 
molecular emission in this region. However, SMA6 was near 
the edge of the ATCA primary beam. 

It is notable that the orientation of the NE-SW outflow in- 
ferred to originate from SMAl is different from the largescale 
SE-NW outflow direction inferred from single-dish SiO(5-4) 
data w ith 23" resolution presented by Megeath & TieftrunkI 
(119991) . In their single dish data, blueshifted emission islo- 
cated to the SE of I(N) in the vicinity of SM2 (Fig. |43] |. 
while blueshifted and redshifted emission overlap in the vicin- 
ity of I(N). Toward I(N), the P-V diagram presented by 
iMegeath & TieftrunkI (119991) shows good agreement with the 
velocity range observed in the SMA SiO(5-4) data (-30 to 
-1-17 km s"', Fig.|7t,b), and indeed the SEST data did not have 
sufficient spatial resoluti on to resolve th e SMA sizescale flow 
from SMAl. Moreover. lMegeath & TieftrunkI (1 19991) also ob- 
serve a clump of shocked 2.2 /j,m H2 emission coincident 
with the northeastern lobe of the SMAl outflow and extended 
4.5 /im emission . Thus, it seems likely that the single outflow 
reported by Meg eath & Ti eftrunk ( 1999), is in fact at least two 
outflows - one consistent with the NE-SW flow detected by 
the SMA and a second in the vicinity of SM2 and the CI and 
C2 water masers. This second outflow would also explain the 
widespread 44 GHz CH3OH maser emission observed in this 
region. While this seems the most probable interpretation of 
the data, it is also possible that spatial filtering of the SMA 
data on sizescales > 20" have played a role in the apparent 
discrepancy. Unfortunately, the SM2 area is beyond the pri- 
mary beam of the SMA observations, so additional data will 
be required to explore this region of activity further 

4.2. Nature of Mid-infrared Emission 

To assess the nature of the mid-IR emission in more detail 
we analyzed the 1-8 /im flux densities of several knots in the 
vicinity of I(N) that morphologicaUy appear to be dominated 
by shock line emission. We included data from th e Spitzer 
IRAC bands (see for example I Ybarra & Ladall2009l) . b ut also 
near-i nfrared J, H, and Ks fluxes (see for example iSmithI 
119951) . To derive the physical conditions, we have imple- 
mente d a set of shock models sirrdlar to that of'Ybarra & Ladaj 
(I2OO9I) (see alsolNeufeld & Yuanl200 8; Smith & Rosen 200ll 
ISmith et al]l2006l) ~ We calculate the equilibrium excitation 
of molecular hydrogen (H2) in a grid of temperatures (1000- 
5000 K) and densities (10 ^-10^ cm~^) using the escape prob- 
ability code Radex (van d er Tak et al. l HOOTl since the optical 
depths are low, the escape probability formalism is not crit- 
ical, but the code was still very useful to read the extensive 
tables of molecular data). W e use excitation rates from col- 
Usions with H2 and He from iLe Bourlot etall (119991 12002!) 
including reactive collisions as prescribed therein, and from 
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Fig. 1 1 . — Spitzer 3-color image with RGB mapped to 24, 4.5, and 3.6 fim. Wliite 1.3 mm continuum and cyan SiO integrated intensity contours are overlaid, 
along with yellow + symbols showing the locations of the 44 GHz CH3OH masers. The 1.3 mm continuum contour levels are the same as Fig. 1, and the SiO 
(5^) integrated emission contour levels are 6, 10, 14, 20, 26 mjy beam"' *km s"'. 



collisions with H from iWrathmall et al.' ('2007'). Quadrupole 
transition probabilities are from Wolniewicz et al. ( 1998). We 
convolve the line intensities with Spitzer (IReach et al.ll2005l) 
and 2MASS'^ filter transmission profiles to calculate band- 
averag e fluxes as specified in appendix A of iRobitaille et al.l 
(l200l . 

The analysis of the shocked emission is complicated by dif- 
fuse PAH emission, especially in the 5.8 and 8.0 pim images. 
We performed several different background subtractions to at- 
tempt to isolate the H2 emission, with the assumption that 
most of the emission at 4.5 fim is H2 (there are no PAH emis- 
sion features in that band). Naturally this decomposition is 
subject to uncertainty, so the precision of shock physical pa- 
rameters is limted by this observational fact. We find that 
the near and mid-lR broadband flux ratios are most consis- 
tent with the hottest post-shock temperatures ('-~^4000K), in- 
dicating that these are strong shocks as would be expected 
from a powerful massive stellar outflow. The density in the 
shock cooling zone is harder to constrain because density 
changes the broadband flux ratios less than does tempera- 
ture (see Fig. 2 in Ybarra & Lada 2009), and because fore- 
ground extinction can mimic a decrease in emitting material 
density. The fluxes here are broadly consistent with densi- 
ties of 10^-10^ cm ^. Once we have constrained the phys- 
ical parameters in the H2 emitting zone, we can translate 
mid-infrared fluxes into total shock luminosity. For shock 
recombination zones (T~4000K) our models indicate that 

'^ |http : / /www, ipac ■ caltech . edu/2mass /releases /allsky/ d 



6% of the total H2 line emission comes out in the IRAC 
4.5 fim band. The NGC63341(N) outflow emits 15±3 Jy over 
100±20 square arcseconds (sum of both lobes). At a distance 
of 1 .7 kpc, that corresponds to 2300±300 Lq total luminos- 
ity in shocked molecular hydrogen, with an average column 
density of 3.5±0.5 x lO^'cm-^, for a total of 0.1 Mq of hot 
shocked material (that which is emitting at the highest tem- 
peratures - a much larger mass of cooler and entrained ma- 
terial is present in the outflow as well, simply not emitting in 
the near infrared). 

SMA4 is the only source in 1(N) associated with 6.7 GHz 
CH3OH masers (see §3.3.2), and in §3.6 we report the detec- 
tion of a 24 fim source 2" to the SW of the 1.3 mm contin- 
uum peak. Class 11 6.7 GHz CH3OH masers are thought to 
be radiatively pumped by warm dust, and are f ound almost 
exclu s ively i n regions of massiv e star formation ( Cragg et alj 
120021 120051; iMinier et al.l l2003h . providing strong evidence 
that SMA4 is a massive protostar. Two scenarios exist for 
the nature of the 24 /im emission: (1) it may trace emission 
from hot dust near a central protostar and is the 24 /im coun- 
terpart to SMA4 or (2) it traces emission from hot dust in 
the walls of the cavity formed by the outflow traced by ex- 
tended 4.5 jum and molecular line emission, as well as 44 
GHz methanol r nasers W/SW of SMA4 . Recent work by 
iDe Buizer et"al] (120051) : iDe Buizeri (12007 ') favor the latter in- 
terpretation. Higher-resolution mid-infrared data would help 
to distinguish between these possibilities for SMA4. The 
__24 /im flux is 0.45±0.04 Jy, or 2 Lq (just in the 24 fim band). 
o i^fe esstajt jp feajrit Tfe u ghly 5% of their luminosity in the 24 /im 
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band over a relatively large range of evolutionary state, so if 
this source represents a significant fraction of the luminosity 
of the driving protostar, the latter would have an approximate 
luminosity of only about 40 Lq. This low luminosity lends 
support to the interpretation that the 24 /im emission is com- 
ing from hot dust away from the central source, such as on an 
outflow cavity, and that the driving source itself is still highly 
embedded and extinguished even at 24 /im. 

4.3. What are the 44 GHz Methanol Masers Tracing? 

The 44 GHz Class I CH3OH maser line is oft en found in 
regions of high-m a ss star formation j Haschick et al.l 119901: 
Slvsh et al.' 19941; iKurtzeLaLl l2004t iVaFTts & Larionovl 
2007; Pratapetal. 2008|~~ Like other Class I masers 
dPlambeck & Mentenlll99 0), it is often found separated from 
compact Hll regions and in many cases appears to be 
associ ated with shocked gas created by b ipolar outflows 
(e.g. Kurtz et al.l l2004t ISandeU et al.l I200I . While fliis 
maser has been found within other millimeter protoclusters 
such as S255N an d G31.41h-0.31 (Cvgano wski etaH l2007t 
lArava et"aLl |2008l ), the number and spatial distribution of 
maser spots is particularly large in NGC 6334 I(N). In all of 
these studies including the current one, the 44 GHz CH3OH 
masers tend to lie within a few km s~' of the systemic velocity 
even though they typically trace powerful outflows. This di- 
chotomy is likely due to the fact that an enhanced column of 
CH3OH gas with sufficient velocity coherence along the line- 
of-sight is most likely to be found where a shock is moving 
perpendicular to the line-of-sight (i.e. a shock viewed edge- 
on). The same effect is believed to give rise to the narrow NH3 
(3,3) maser emission located near the ends of the protostellar 
jet from IRAS 20126H-4104 (Zhang et al. 1999). 

As described in §4.1, it is possible that the 44 GHz masers 
located in the SE portion of the observed field of view are 
associated with an outflow originating from the relatively un- 
explored region of continuum emission SM2. The discov- 
ery of an X-ray source near SM2 indicates that additional 
cluster members may be present in this area. A more gen- 
eral conclusion can be drawn from the relatively widespread 
extent (45" 1.2 x 10'^ cm) of 44 GHz maser emission, 
an extent that is shared by ammonia (1,1) and (2,2) line 
and 450 /m continuum emission (see Fig. ). For a sin- 
gle dominant YSO to maintain this region at a conserva- 
tive minimum t empera t ure of 30 K (the ammonia tempera- 
ture from Kuiper et al.l (Il995h and average dust temperature 
from ISandelU (120001) 7 would require a central luminosity of 
5 X 10^ Lq (e.g. iScoville & Kwanlll976 ) which is 30 times 
higher than the value (1.7 x lO'* L©) inferred from the far- 
infrared spectral energy distribution for this region (Sandell 
12000) . Thus, these phenomena must be driven by the collec- 
tive activity of a number of cluster members. 

4.4. Temperatures of the Hot Core Sources SMAl and SMA2 

When several transitions of the same molecule spanning 
a range of energies are observed, temperatures can be in- 
ferred from the rotation d iagram method (see for example 
[Goldsmith & Langeill999h . using the relations 

— = ^-^ 2^ / S^dv, (1) 

and 

logiNjgu) = logiN,o,/Q(Tro,))-0A343E,/mo„ (2) 
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Fig. 12. — Rotation diagrams for tlie CH3OH transitions observed for 
SMAl (top) and SMA2 (bottom); the fitted temperatures are indicated. The 
column densities have been corrected for optical depth effects as described in 
§4.3. 

where k is Boltzmann's constant, ly is the frequency, gj and 
gK are the degeneracies associated with the nuclear spin and 
k quantum number respectively, fi^ is the square of the dipole 
matrix element, S is the statistical line strength, / Si^dv is the 
observed integrated intensity of the transition. A',,,, is the total 
column density, Q{T,-oi) is the partition function evaluated at 
the rotation temperature Trot, and £„ is the upper state energy 
of the transition. There are seven transitions of CH3OH de- 
tected toward SMAl and six toward SMA2 (including both A 
and E type transitions) that can be used for this purpose. The 
/ S,jdv for each transition was calculated from Gaussian fits 
to the line emission at the 1 .3 mm continuum peaks. The rota- 
tion temperature resulting from a weighted least squares fit to 
the data is 189 ± 19 K for SMAl and 154 ib 14 K for SMA2 

However, as demonstrated by Golds mith & LangeJ ( 11999b 
and others, line optical depth can artificially increase the de- 
rived rotation temperature if not corrected for. In general the 
line optical depth in LTE can be calculated from 

r = ^^^^,^Se-^J^He'-l"- - 1)-^, (3) 
3h go Av 

where h is Planck's constant. No/ go is the column density of 
the ground state divided by the ground state degeneracy, T^x is 
the excitation temperature, and Av is the FWHM line width. 
Then the left hand side of Eq. 2 is modified to log{CrNu/gu) 
where Cr = t/(1 -e~^). Since we do not know the ground 
state column density, this equation cannot be used directly 
to calculate r for each transition. Unfortunately, in the I(N) 
SMA dataset we do not have a suitable isotopologue that can 
be use d to infer the CH3O H optical depth either (see for ex- 
ample Bro gan et al.l2007b . Instead, we have iterative ly solved 
for the optical depth and Trot that produces the best fit to 
the data (i.e. minimizes the scatter). Using this technique 
we find the best fit occurs for a Troi=l65±9 K for SMAl 
and Tro,=l45 ± 12 K for SMA2 (see Figure [Eli), and mod- 
est optical depths of 1.75 for SMAl and 0.5 for SMA2 for 
the observed transition that would have the highest optical 
depth (CH3OH-E (4+2,2-3+1,2)) at 218.4400 GHz). While 
not dramatically different from the Trot derived without opti- 
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cal depth correction, these corrected temperatures are clearly 
more accurate. We do not have enough transitions to do a 
more sophisticated non-LTE analysis or to compare the re- 
sults for A and £ transitions independently (see for example 
iSutton et al']|2004 >). However, we note that for the CH3OH-E 
(5+1 4-4+2,2) tran sition (at 216 . 9456 GHz) which we have in 
common with the lSutton et al.l (120041) methanol survey of the 
hot core source W3(OH) with comparable spatial resolution, 
the derived opacity is similar for the two sources (0.45 for 
SMAl in 1(N) compared to 0.65 for W3(OH)). Indeed the de- 
rived temperatures for SMAl and SM A2 are also very similar 
to that of W3(OH) (r,o, = 140 ± 10 K, ISutton et al . 2004). 

The optical depth implied by the brightness temperature 
of the strongest methanol transition (4+2,2 - 3+i, 2, 24 K) to- 
ward SMAl is only 0.16 assuming that 7i,„=rrof=165 K and 
the emitting region is the size of the LSB synthesized beam 
(2"39 X 1"55). An optical depth of 1.75 for this transition 
implies that the emitting region is actually only 0"8 in size 
(1400 AU). This size is in good agreement with the fitted size 
(0"9 X 0."7, PA=162°) of the integrated intensity measured 
from the compact, but still strong CH30H-A''"(16i,i6- 152, 13) 
ti-ansition (£/ = 219.8 cm-')at 227.8147 GHz. This size im- 
plies that the hot core region encompasses both SMAlb and 
SMAld(seeFig.|2]3). 

The column densities of CH3OH for SMAl and SMA2 are 
(1.0±0.4) X 10'^ cm-2 and (5.6 ±0. 6) x lO"^ cm"^ re spec- 
tively, an order of magnitude less than lSutton et al.l (120041) find 
for W3(OH). Sutton et al. ( 2004) find a CH3OH abundance of 
'--^ 2 x lO"*" for W3(OH). Usi ng single dish obse rvations of 
infrared dark clouds (IRDCs). Leurini et al.l (I2007h find abun- 
dances an order of magnitude smaller for the IRDC cores. Us- 
ing this range of abundance and assuming that the hot core 
size is 1400 AU, we estimate that the H2 column density is 
(0.5-5) x 10^3 cm-- and the H2 density is (0.2-2) x 10^ cm'^ 
for SMAl and about a fac tor of two sma lle r for S MA2. 

Using the ATCA, iBeuther et all (120081) observed 
CH3CN(5if-4^) up to K=4 toward SMAl and SMA2 
and derived a rotation temperature of 170 ±50 K for 
SMA2 in good agreement with our CH+30H measurement 
(145 ± 12 K). The optical depth in SMAl was too high to 
obtain a CH3CN Trot for SMAl in good agreement with our 
finding that this source has a hig her column den sity than 
SMA2. Also using the ATCA, Beu ther et al] (120071) detected 
NH3 (5,5), and (6,6) emission toward both SMAl and SMA2 
confirming that is greater than 100 K, though it was not 
possible to derive more accurate estimates from those data. 

4.5. The Nature of the 1.3 mm Sources 
4.5.1. Spectral Energy Distributions 

Because no higher frequency arcsecond-resolution data ex- 
ist on this region, the determination of the nature of the 

1.3 mm sources hinges entirely upon accurate measurements 
of their spectral index toward longer wavelengths (thus mo- 
tivating our careful re-reduction of the available data). The 
comparatively poor angular resolution of the current 1.3 and 

3.4 mm data compared to the 7 mm, 1.3 cm, and 3.6 cm data 
makes it challenging to disentangle the free-free versus dust 
contributions to these sources. Figure [13] shows the spec- 
tral energy distributions (SEDs) of SMAl, SMA4, and SMA6 
based on the flux densities provided in Table [T] The different 
resolutions of the different datasets are definitely a concern. 
However, in many cases the emission does not appear to be 
resolved even in the highest resolution observations of SMA4, 



and to some degree SMA6. In that case, as long as no signifi- 
cant emission is resolved out, the data will still be comparable 
despite the differing resolutions. Caveats for a few sources are 
described in detail below. 

For SMAl, which contains at least four components at the 
longer wavelengths, we show both the integrated emission 
over the whole SMAl region, and SMAlb by itself at wave- 
lengths longer than 7 mm for comparison. The integrated 
emission from SMAl shows evidence for both a free-free and 
an optically thin dust component. The integrated dust spec- 
tral index of adust = 3.8 ±0.1 {Su oc v"') was determined by a 
linear least-squares fit to the three flux densities from 7 mm 
to 1.3 mm (for SMAl at 7 mm this is the sum of emission 
from SMAla,b,c, and d; the free-free contribution is an or- 
der of magnitude smaller and hence negligible in the dust fit). 
On smaller scales, we also find that the 7 mm emission from 
SMAlb is dominated by dust, while the 3.6 cm emission is 
consistent with optically-thin free-free emission with a spec- 
tral index of i^"*", such as that from an Hll region with an 
electron density of 3 x lO'* cm~^, temperature of 10000 K and 
a diameter of 0.3" (510 AU). At the intermediate wavelength 
of 1 .3 cm, these two emission mechanisms are comparable for 
SMAlb, with « 30% of the emission arising from dust. This 
is a notable result since 1.3 cm data are typically assumed to 
be free of dust. 

The emission from SMAl a may all be due to dust since it 
is not detected longward of 7 mm; the 3cr 1.3 cm upper limit 
of ^ 0.18 mJy implies that the spectral index must be > 2.4. 
In contrast, both SMAlc and SMAld are detected at 1.3 cm 
and have a 1.3 cm to 7 mm spectral index of siml. This spec- 
tral index could be due to optically thick dust emission which 
seems very unlikely at these wavelengths, optically thick free- 
free emission such as one might observe around a very dense 
hypercompact H II region, or a combination of optically thin 
free-free and dust emission. If we take the total observed flux 
density of SMAl at 1.3 cm (1.07 mJy, i.e. the sum of SMAl a 
upper limit and SMAlb,c,and d detections) and subtract the 
total dust model (0.50 mJy) and the SMAlb free-free emis- 
sion model (0.35 mJy), we find that 0.22 mJy remains, sug- 
gesting that the integrated dust model cannot account for all of 
the SMAlcH-SMAld emission. Thus, one or both of SMAlc 
and SMAld must have a weak free-free component, and since 
these objects have 1.3 cm flux densities of 0.27 and 0.30 mJy, 
respectively (greater than the residual), both must also have 
a dust component. If all of the free-free residual (0.22 mJy) 
belonged to either SMAlc or SMAld, the 3cr 3.6 cm upper 
limit of ~ 0.12 mJy implies a free-free spectral index > 0.6, 
consistent with a thermal jet, stellar wind, or hypercompact 
H II region. If instead the residual free-free emission is about 
equally split between SMAlc and SMAld both could harbor 
small optically thin H II regions in addition to a dust compo- 
nent. 

Understanding the nature of SMA4 is complicated by the 
lack of positional agreement between the 1.3 - 3.4 mm data 
compared to the longer wavelength data (3.6 cm, 1.3 cm, and 
7 mm offset ^ 0."4 to the NW; Fig. |3]l. From the spectral 
break in the SMA4 SED around 7 mm shown in Fig. [T3]it is 
clear that both dust and free-free emission is present in this 
region. The spectral index of the SMA4 emission between 
1.3 and 3.4 mm is adust = 4.4 ±0.3, while the spectral in- 
dex between 3.6 cm and 7 mm is a j-j- = 0.63 ± 0.25 suggest- 
ing either a hypercompact H II region wi th a turnover a round 
7 mm (see for example the review by lLizanoll2008l) . or a 
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Fig. 13. — Spectral energy distribution.s for sources SMAl, SMA4. and SMA6. For SMAl, the sum of the emission of SMAla.b.c, and d are shown as solid 
squares, while the flux density as a function of frequency for SMA lb is shown as open squares for u <43 GHz. Open squares are also used for the u <43 GHz 
emission for SMA4 to indicate that they may not be from the same source as the higher frequency data, the fits in this case assume case 3 described in S4.5.1. 
Emission models from dust are shown as dashed (or dot-dash in the case of SMAlb) lines, while emission models for free-free emission are shown as dotted 
lines. For SMAl the sum of the dust and free-free emission models, described in detail in §4.4.1, is shown as a solid line. 



thermal jet/stellar wind ( lRevnoldslll986l) . Overall, the data 
are consistent with the following three possibilities: (1) the 
emission all arises from the same source which contains both 
dust and a hypercompact H n region (i.e. the position off- 
sets are just within the combined position uncertainties); (2) 
Two sources are present: one is a compact dust source with 
no detectable free-free emission (1.3 - 3.4 mm emission) and 
the other (7 mm - 3.6 cm emission) is a hypercompact Hn 
region (> 700 AU away); or (3) the free-free emission arises 
from a one-sided thermal jet driven by the dust source of case 
(2). If we assume case (1) and fit the centimeter wavelength 
emission from SMA4 after extrapolating and subtracting the 
dust spectrum, we find that a free-free component with a con- 
stant density, temperature, and size of 1 .2 x 10^ cm"-', 8000 K, 
and 53 AU best approx imates the 3.6 and 1.3 cm data points 
(see iHunter et alj|2008l for the details of the free-free emis- 
sion model). However, this model overpredicts the 7 mm 
flux density by a significant margin. More sophisticated free- 
free emission models, such as those with a density gradient 
and/or a shell geometry, may provide a bette r match to the 
7 mm data (see for example lAvalos et"al]|2006l) . However, the 
present data are insufficient in angular resolution and number 
of spectral measurements to meaningfully constrain a more 
detailed model. For case (2), it is even more difficult to fit the 
7 mm data with a hypercompact H 11 region model because the 
spectrum from 3.6 cm to 7 mm is so linear, and yet not steep 



enough to still be optically thick (i.e. spectral index of 2). 
Thus, case (3) seems to best fit the current data (see Fig. [T3] 
with the 3.6 cm to 7 mm data having a typical jet or stellar 
wind spectral index (ajj- = 0.6 Revnolds 1986), and no dust 
component. We note that in the jet interpretation for the cm- A 
emission (case 3) the ~ 0."4 offset of the cm-A to the west is in 
good agreement with the fact that a westward outflow appears 
to emanate from SMA4. 

Although there is an offset ^ 0"3 to the north between the 
7 mm and 1.3 mm emission toward SMA6, the fact that the 
shape of the continuum contours at these two wavelengths is 
similar, and most of the apparent offset is along the long axis 
of the 1.3 mm beam makes it quite plausible that they are in 
fact coincident (also note that the offset is not in the same di- 
rection as that of SMA4). Indeed, the emission from SMA6 
can be explained completely by optically-thin dust emission 
(at the current level of sensitivity). The dust spectral index 
from a linear fit to the three flux densities from 7 mm to 
1.3 mm is a = 4.0 ±0.2. If we assume all the 1.3 mm emission 
arises from a region with the size measured from the higher 
resolution 7 mm image, (0"9 x 0."7), the brightness tempera- 
ture is Th = 25 K setting a lower limit on the dust temperature. 
There is a hint of a break towards a flatter spectral index at 
1.3 mm for SMA6, as would be expected if the dust emission 
were becoming optically thick, but higher frequency data are 
required to confirm this trend. 
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SMA2 and SMA3 are the only strong, compact millime- 
ter sources not to be detected at wavelengths longward of 
3.4 mm. Between 1.3 and 3.4 mm the spectral index of SMA2 
is a = 4.0 ± 0.3, and that of SMA3 is a = 4.3 ± 0.3, consis- 
tent with optically thin dust emission. Extrapolating these 
spectral indices to 7 mm, we would have expected to detect 
SMA2 and SMA3 at the 1.4 and 0.9 mJy level, respectively 
(amounting to la and 4.5ct) if the emission were unresolved 
at 7 mm. The sizes derived for SMA2 and SMA3 at 1.3 mm 
are 2"1 X 1 "3 and 2. "5 x 1 "9. If these are the true sizes of the 
emitting regions, the peak 7 mm intensity would be dimin- 
ished by a factor of ~ 10 over that predicted here, plausibly 
explaining the non-detections. The emission denoted VLA3 
is only detected at 7 mm, and thus its nature is very uncertain. 
The non-detection of VLA3 at 1.3 cm implies that its spectral 
index is > 2, and we estimate a 1.3 mm upper limit of 150 
mJy suggesting that its spectral index must also be less than 
~ 3. These limits suggest that VLA3 is predominately due 
to dust emission, but contribution from a hypercompact H II 
region cannot be excluded. 

In summary, of the nine compact sources detected in I(N), 
all appear to have a dust component. Four of these sources 
(SMAlb, SMAlc, SMAld, and SMA4) also show evidence 
for free-free emission at our current level of sensitivity. Three 
sources (SMAlb, SMAld, and SMA2) and possibly a fourth 
(SMA4) show modest hot core molecular line emission, also 
suggesting the presence of a central powering source. One 
source with neither free-free nor hot core emission does ap- 
pear to be powering an outflow (SMA6). From this circum- 
stantial evidence it appears that at least six of the compact 
so urces may harbor a central source. However, as described 
by IZinnecker & York3 (120071) . in deeply embedded regions 
it is extraordinarily difficult to determine whether a massive 
protostar or zero age main sequence star (ZAMS) has formed, 
as this designation is dependent upon detecting the source in 
the near-lR (a wavelength regime in which 1(N) is still ob- 
scured). In principle, the detection of free-free emission is 
often taken as indicative of at least an early B-type star How- 
ever, even the presence of free-free emission is ambiguous 
since the accretion luminosity itself can be sufficient to ionize 
the gas around a protostar. Thus, while it is clear that a num- 
ber of massive to intermediate stars are forming in 1(N), the 
exact number and evolutionary state is uncertain. 

A few comments on the previously published flux measure- 
ments for these sources are in order When the more ac- 
curate absolute flux calibration presented here for the orig- 
inal ICarral et al. (2002) 3.6 cm data is taken into account, 
the agreement for SMA4 is quite good. However, it is un- 
clear why the detection of SMAlb was not reported by these 
authors, since it is 1.6 times brighter than SMA4. We ob- 
tain 7 mm flux densities between 2.2 (SMA6) to 8 (SMAla) 
times smaller than those reported by Rodriguez et al. ( 2007), 
though our 1.3 cm flux densities are comparable. This dis- 
crepancy comes in spite of the fact that the flux densities 
derived for the phase calibrators at both 7 and 1.3 cm are 
within 2% for o ur red uction compared to that reported in 
iRodrfguez et al.l (l2007h . As described in §2.2.3, continuum 
emission from the strong ultracompact H II region in source 
1 causes significant contamination of the 1(N) region if not 
included in the cleaning; however our detailed tests suggest 
that at most it could erroneously increase the integrated fluxes 
by a factor of ~ 2, and thus cannot account for the large dis- 
crepancy. In addition, the integrated flux densities reported 
in .Rodriguez et al.. (.2007.) do not appear consistent with the 



contour levels plotted in their Figures 2 and 3, indeed the 
peak values are reasonably consistent with those reported here 
so the problem may have occurred in their integrated Gaus- 
sian fits. Finally, o ur re-imaging of the 3.4 mm data from 
iBeuther et al.l (120081) has increased the peak flux densities by 
a factor of about 1.5 due to our self-calibration and correction 
for the primary beam attenuation. 

4.5.2. Mass estimates from the dust emission 

In Table|6]we give estimates for the gas mass Mgas, molec- 
ular hydrogen column density Nh2, and molecular hydrogen 
density mh, for all seven 1 .3 mm dust sources. It should be 
noted that if a protostar or zero age main sequence star has 
already formed, the millimeter data is only sensitive to the 
circum-(proto)stellar material. The gas masses were calcu- 
lated from 



3.24 X 10-35i,(Jy)Z)2(kpc)iJC^„„, 



Mr, 



(4) 



Jiy, 7j„s,)i/3(THz)Ky 

where Si, is the integrated flux density from Table [1] D is the 
distance (assumed to be 1 .7 kpc), R is the gas to dust ratio (as- 
sumed to be 100), J{.v,T) = i/iexpihiy/Tdustk)- 1), and Q,„/ 
is the correction factor for the dust opacity Crj„„ = Tdust/i^ - 
e"^''""). We have used mm = 1 c m^ g~' which is the av- 
erage of the dust opacities derived bv lOssenkopf & HenningI 
(.1994) for the case of thin ice mantles and densities between 
10*^ to 10^ cm-3 at 1.3 mm. We have estimated a range of dust 
temperatures based on the observed spectral line properties 
of each source: we use the derived Trot values calculated in 
§4.4 for SMAl and SMA2, and a lower temperature range for 
SMA3 to SMA7 based on the non-detection of hot core spec- 
tral lines towards these sources. The column densities were 
estimated from 



2.4 X io"'X(Jy)/^c^,„,„ 
J{y,Tdus,)vH^'^7)0i(")02O^iu 



cm 



(5) 



where 9\ and 6*2 are the fitted sizes from Table[T] The nn^ was 
calculated by dividing Nh^ by the geometric mean of the fit- 
ted size, assuming a distance of 1.7 kpc. The estimat e d dus t 
masses are similar to those reported in Hun ter et al.l (120061) 
when the differences in absolute flux calibration and assumed 
dust temperatures are taken into account. As with all such es- 
timates, the Mgas, Nh^, nfj, reported in Tableware uncertain 
to at least a factor of 2. The total mass estimated for the com- 
pact emission of I(N) is between 89 and 232 Mq. These dust 
emission based masses do not include the mass of any pro- 
tostars or ZAMSs that may already have formed within the 
compact cores. It is also notable that [Matthews et al] ( l2008h 
estimate that the mass of the whole I(N) core is ~ 2780 Mq, 
so there is still a considerable reservoir of material that has 
not yet collapsed. 

The values for SMAl are particularly uncertain because 
the current 1.3 mm resolution is insufficient to distinguish 
between the four dust sources detected at 7 mm (SMAla, 
SMAlb, SMAlc, and SMAld). This unresolved structure 
makes the apparent size of this source quite large - signifi- 
cantly affecting all derived quantities. For example, if we as- 
sume that the majority of dust emission comes from a region 
half the size (in both directions) of the current fit (2 ."8 x 1 "9), 
the mass increases only a little (due to increase in Cr^,,,,), 
but Nh-, and both increase by about an order of magni- 
tude. An additional source of uncertainty is the Tdust- As de- 
scribed in §4.3, the hot core line emission appears to arise 
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from a 0"8 region encompassing only SMAlb and SMAld. 
Thus, the Tdust range derived from the CH3OH rotation di- 
agram may not be appropriate for SMAla or SMAlc. Us- 
ing K7 mm = 0.05 cm^ g"' (obtained by scaling K1.3 = 1 by 
(7 mm/1.3 mm)"'^ with (3 = ajust-'^- where a^ust = 3.8 from 
Fig. \T3[. and Tdusi = 155- 175 K we find that the sum of the 
masses derived for SMAla,b,c, and d at 7 mm is 20 - 18 Mq 
in good agreement with the 1.3 mm estimate for SMAl (16 
- 14 Mq). Indeed, the agreement is nearly perfect if we as- 
sume that the dust emission at 1.3 mm arises from an area 5 x 
smaller than the current size estimate (i.e. Table [T]!, which 
is certainly plausible. Higher resolution 1.3 mm observations 
are required to better constrain the temperatures of SMAla 
and SMAlc, as well as the sizes of all four SMAl sources. 

4.6. Evidence for Accretion and Infall Around SMAl ? 

As described in §3.4, a SE-NW velocity gradient of 
5 km s"' is observed toward SMAl in some (but not all) 
of the species with compact emission (a similar gradient 
was inferred from a double peaked NH3(6,6) profile by 
iBeuther et alll2007l) . Since the gradient is barely spatially 
resolved with the current SMA resolution, the nature of this 
gradient is unclear Its direction is more or less perpendic- 
ular to the large scale NE-SW outflow suggestive of a pos- 
sible disk interpretation, but it also has the same orienta- 
tion as a line joining the SMAlb and SMAld sources, and 
thus could simply represent unresolved emission from these 
two sources if they have different systemic velo cities (as was 
found for CephA-East in a similar situation by iBrogan et alJ 
l2007h . We note that this velocity gradient was not observed 
in the CH3CN(5a:^a') transitions observed by IBeuther et al.l 
(|2008) with the ATCA. Higher spatial resolution data will be 
needed to further explore this velocity gradient. 

As described in §3.5.2, there is a trend of increasingly red- 
shifted absorption with increasing line excitation temperature 
toward SMAl, using CN (£,=3.8 cm"'), H2CO(£/=7 cm"'), 
and H2CO (Ei=40 cm"') as tracers (see Fig.[TOli). In §4.3 we 
find that the hot core of SMAl has a temperature of 165 ±9 K, 
demonstrating that this source is definitely centrally heated, 
suggesting higher excitation lines will be found closer to the 
center We suggest that the observed redshift/temperature 
trend is tracing infall that is accelerating with depth into the 
core. It would be very interesting to follow-up this tentative 
result with higher angular and spectral resolution. In partic- 
ular if the continuum emission were resolved, one could be 
assured that all emission/absorption originates from only the 
front side of the source, removing any ambiguities (see for 
example iChandler et a l. 2005). 

4.7. The Location of the Protocluster 

Finally, it is interesting to note that the massive protocluster 
appears to be forming near the northwestern edge of the par- 
ent cloud rather than toward the center. This offset is ap parent 
in both the JCMT 450 /im continuum image (Fig. 14.31 1 and in 
the ammonia (1,1) and (2,2) images tBeuther et al.ii2005i) . It 
is notable that the archival SCUBA 450 im\ images recently 
produced by IDi Francesco et akl ((2008') show very good po- 
sitional agreement (tojwithin 3") with the earlier images pro- 
duced by Sandell (2000) (and shown in Fig. 14. 3b . The ridge of 
dust emission between the protocluster and the masers asso- 
ciated with the SM2 continuum source coincides with a dark 
area in the Spitzer 3-color image (Fig. [TTT i. i.e. an infrared 
dark cloud (IRDC; see for example .Egan et al.,.1998i) . Clearly 



this region is optically-thick in the mid-infrared, and contains 
a significant amount of potential star-forming material. Future 
observations with broader bandwidth and greater sensitivity 
with the EVLA and ALMA may reveal whether this area is 
truly quiescent or contains embedded, pre-protostellar cores. 

5. CONCLUSIONS 

Our multiwavelength radio through near-infrared investiga- 
tion of the massive protocluster NGC6334I(N) provides new 
morphological and quantitative detail on the star formation 
activity of the individual massive cores. We show that I(N) is 
undergoing copious massive star formation with many of the 
dominant sources detected only at millimeter wavelengths. 
The principal source SMAl is resolved into four components 
of dust and/or free-free emission (SMAla, b, c, and d). De- 
lineated by centimeter continuum and water masers, SMAlb 
and SMAld form a close central pair (> 800 AU) which 
powers a hot core (7'ro,= 165 ±9 K) of molecular line emis- 
sion and likely drives the dominant large-scale bipolar out- 
flow which exhibits evidence of precession. We find a spa- 
tially unresolved ~ 5 km s"' velocity gradient across SMAl 
(encompassing all 4 subarcsecond components). The origin of 
this velocity gradient (rotation, or cluster kinematics) requires 
higher resolution millimeter wavelength follow-up. We also 
detect possible infall that is accelerating with depth into the 
SMAl core; higher angular resolution is also required to con- 
firm this finding. The mass of gas around SMAl as traced by 
compact emission is estimated to be ^ 15 Mq, but if this core 
is still actively accreting, this may grow significantly larger in 
the future. 

A weaker hot core (7;o,=145 ± 12 K) is found toward SMA2 
which also exhibits redshifted CN absorption along with with 
highly blueshifted water masers, suggesting simultaneous in- 
fall and outflow in this source. SMA4 consists of a dust core 
with a probable centimeter jet, water maser emission, and 
blueshifted SiO emanating to the west suggesting an outflow 
which terminates in a reflection nebula traced by the 24 //m 
source. Mostly lacking in molecular lines, SMA6 exhibits 
a centimeter-millimeter SED consistent with pure dust emis- 
sion and drives an outflow traced by SiO, CO and 44 GHz 
methanol masers. Overall, the narrow velocity range and spa- 
tial arrangement of the 44 GHz methanol masers with respect 
to the protocluster members are consistent with their origin in 
low-velocity gas impacted by outflows from these objects. 

The presence of multiple outflows in the region suggests 
that many protostars are forming simultaneously. At the same 
time, the variety of phenomena associated with the individual 
cluster members suggests a spread in either age or total lu- 
minosity. Also of interest is the spread in systemic velocities 
observed between the cluster members, this spread is at least 
2.5 km s~' and may provide an important clue about the dy- 
namics of cluster formation. A more quantitative assessment 
of these key properties will require broader wavelength cov- 
erage to better constrain the SEDs, and more sensitive contin- 
uum and line observations enabled by the EVLA and ALMA 
in the near future. At the same time, a deeper exploration 
of the mostly featureless, infrared-dark region between the 
protocluster and SM2 may prove fruitful in the study of pre- 
protostellar cores. 

We thank Goran Sandell for providing us with the calibrated 
JCMT 450 /im image. This research has made use of the 
NASA's Astrophysics Data System Bibliographic Services, 
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TABLE 1 

Observed Properties of Compact Continuum Sources 





Centroid position" 


Peak 


Integrated 


Estimated 


Source 


RA (J2000) 


Dec (J2000) 


Intensity 


Flux Density 


Size 




Q,(h m 


5{° ' ") 


(mjy beam"' ) 


(mJy) 


(" X " [P.A. °]) 








1.3 Jvtillinieter^ 






SMAl 


17 20 55.188 


-35 45 03.92 


826 (6) 


2550(50) 


2.8 X 1.9 [23] 


SMA2 


1 / 20 54.890 


-35 45 06. / I 


324 (6) 


980 (25) 


2.1 X 1.3 [137] 


SMA3 


17 20 54.990 


-35 45 07.29 


258 (6) 


1060(31) 


2.5 X 1.9 [9] 


SMA4 


17 20 54.667 


-35 45 08.38 


141 (6) 


770 (40) 


3.0 X 1.5 [41] 


SMA5 


17 20 55.054 


-35 45 01.49 


140 (6) 


440 (100) 


< 2.0 


SMA6 


17 20 54.563 


-35 45 17.71 


282 (9) 


600 (29) 


2.2 X 0.9 [38] 


SMA7 


1 7 90 Szl Q^l 
i / ZU V J I 


ac A A S7 7Q 


60(6) 


360 (40) 


3.2 X 2.3 [47] 








3.4 MtUimeter"'^ 






SMAl 




ac Az 7Q 


55 (5) 


90 (10) 


<3.0 


SMA2 


17 20 54.86 


-35 45 06.5 


24(5) 




< 1.6 


SMA3 


1 7 90 S4 QQ 


AS 07 


20(5) 




< 1.6 


SMA4'' 


1 7 on '\A fxCi 

i / zU j4.oy 


ylC A7 C 


13(3) 




< 1.6 


SMA5'' 


17 20 55.04 


-35 45 01.3 


15 (3) 




< 1.6 


SMA6 


17 20 54.588 


-35 45 17.86 


31(3) 




< 1.6 








7 Millimeter 






SMAla 


1 T OA CC 1 CA 

17 zU 55.150 


OC AC AC oo 

-35 45 U5.zis 


0.9 (0.2) 


1.7 (0.3) 


< 1 


SMAlb 


1 7 OA <C 1 Qyl 
1 / ZU J J . 1 7*+ 


ac yi^ n^ sn 


1.3 (0.2) 


2.5 (0.4) 


0.6 X 0.4 [130] 


SMAlc 


1 1 OA C< OT^ 

1 / zU DD.z/D 


"^C /1C AO Af^ 


0.9 (0.2) 




<0.4 


SMAld'' 


17 20 55.229 


-35 45 03.68 


1.0(0.2) 




<0.4 


VLA3'' 


17 20 54.999 


-35 45 06.48 


1.0(0.2) 


1.8 (0.3) 


< 1 


SMA4'' 


17 20 54.624 


-35 45 08.24 


0.7 (0.2) 




<0.4 


SMA6 


17 20 54.585 


-35 45 16.97 


1.3 (0.2) 


1.7 (0.3) 


0.9 X 0.7 [15] 








1.3 Centimeter 






SMAla 






< 0.18 (0.06) 






SMAlb 


17 20 55.191 


-35 45 03.53 


0.50 (0.07) 




<0.3 


SMAlc 


17 20 55.270 


-35 45 02.45 


0.27 (0.06) 




<0.4 


SMAld'' 


17 20 55.229 


-35 45 03.68 


0.30 (0.06) 




<0.3 


SMA4 


17 20 54.623 


-35 45 08.25 


0.47 (0.07) 




<0.3 


SMA6 






< 0.18 (0.06) 












3.6 Centimeter'' 






SMAla 






< 0.12 (0.04) 






SMAlb"* 


17 20 55.202 


-35 45 03.74 


0.39 (0.04) 




<0.3 


SMAlc 






< 0.12 (0.04) 






SMAld 






< 0.12 (0.04) 






SMA4 


17 20 54.623 


-35 45 08.25 


0.25 (0.04) 




<0.3 


SMA6 






< 0.12 (0.04) 







" Relative uncertainties are indicated by the number of significant figures shown. Absolute uncertainties are at least an order of magnitude larger. 

Because the data have different resolutions, the peak intensity values should only be compared for cases where the source is unresolved (indicated by • • • sym- 
bols in Integrated Flux Density column). Estimated statistical uncertainties indicated in parenthesis. 

Fits based on uniformly weighted image. SMAl, SMA4, and SMA6 appear to also have a more compact core so the fitted peaks are underestimated. 

New detection based on archival data. 

Primary beam correction makes a significant difference for these images; the original papers do not appear to have performed this correction based on reported 
flux densities. 

' Extra uncertainty for the peak and integrated fluxes for SMAl and SMA5, as well as SMA2 and SMA3 were added because these source pairs were barely 
resolved from each other. 
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TABLE 2 

22.235 GHz WATER Maser GROUPS IN NGC 6334 I(N) 





Centroid position" 


Associs-tcd 


Peak 


Peak 


Velocity 


# 


RA(J2000 ) 


Dec (J2000) 


mni source 


Intensity 


Velocity 


Range 




a ™ 


(5(° ' " ) 






(km s"') 


(km s^^ ) 


1 


17 20 56.058 


-35 45 32.59 


SM2'' (4.50") 


38.3 


-6.0 


+5.2 to -16.2 


2 


17 20 55.650 


-35 45 32.48 


SM2'' (1.78") 


9.2 


-9.0 


+7.2 to -15.6 


3 


17 20 54.600 


-35 45 17.27 


SMA6 (0.67") 


8.8 


-2.2 


+10.4 to -12.6 


4 


17 20 54.152 


-35 45 13.66 




29.0 


-2.7 


-0.7 to -4.4 


5 


17 20 54.618 


-35 45 08.65 


SMA4(2.11") 


17.7 


+0.2 


+12.4 to -12.9'= 


6 


17 20 54.870 


-35 45 06.32 


SMA2 (0.55") 


2.7 


-13.3 


-4.7 to -15.6 


7 


17 20 54.819 


-35 45 06.18 


SMA2 (1.11") 


5.3 


-0.7 


+0.2 to -1.7 


g 


17 20 55.192 


-35 45 03.77 


SMAlb (0.24") 


78.1 


-6.3 


+13.1 to -26.8'' 


9 


17 20 55.228 


-35 45 03.66 


SMAld (0.02") 


75.7 


-1-3.9 


+5.2 to +2.2 


10 


17 20 55.147 


-35 45 03.78 


SMAl (0.53") 


3.4 


-8.6 


-7.7 to -9.6 


11 


17 20 55.215 


-35 45 02.02 


SMAl (1.96") 


5.5 


-AA 


-3.7 to -6.3 



" Intensity weighted positions for masers listed under each group name in Table|3] 
'° See Sandell (2000). 

'= Emission extends beyond the positive velocity extent of the bandpass. 
'' Emission extends beyond the negative velocity extent of the bandpass. 
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TABLE 3 

22.235 GHz Water Masers Detected in NGC 6334 1(N) 



Velocity 




Centroid position 




Peak 


(km s-i) 


R.A. (J2000) 




Dec (J2000) 


Intensity 








Si° ' ") 


(Jy beam"' ) 






Group CI 






+5.2 


17 20 56.0583 




-35 45 32.4934 


0.384 


+4.9 


17 20 56.0599 




-35 45 32.5678 


0.548 


+4.5 


17 20 56.0585 




-35 45 32.5920 


0.586 


+4.2 


17 20 56.0578 




-35 45 32.5805 


0.893 


+3.9 


17 20 56.0581 




-35 45 32.5797 


1.373 


+3.5 


17 20 56.0586 




-35 45 32.5776 


1.764 


+3.2 


17 20 56.0579 




-35 45 32.6134 


1.675 


+2.9 


17 20 56.0580 




-35 45 32.5594 


1.534 


+2.5 


17 20 56.0581 




-35 45 32.5946 


1.185 


+2.2 


17 20 56.0584 




-35 45 32.5566 


0.934 


+1.9 


17 20 56.0584 




-35 45 32.5619 


0.639 


+1.6 


17 20 56.0614 




-35 45 32.4557 


0.508 


+1.2 


17 20 56.0617 




-35 45 32.4066 


0.381 


+0.9 


17 20 56.0555 




-35 45 32.5404 


0.535 


+0.6 


17 20 56.0579 




-35 45 32.5826 


0.504 


+0.2 


17 20 56.0590 




-35 45 32.5242 


0.497 


-0.1 


17 20 56.0607 




-35 45 32.4663 


0.422 


-0.4 


17 20 56.0584 




-35 45 32.5610 


0.400 


-0.7 


17 20 56.0597 




-35 45 32.5347 


0.558 


-1.1 


17 20 56.0566 




-35 45 32.6202 


0.652 


-1.4 


17 20 56.0588 




-35 45 32.5413 


0.578 


-1.7 


17 20 56.0581 




-35 45 32.6271 


0.749 


-2.1 


17 20 56.0587 




-35 45 32.5695 


1.507 


-2.4 


17 20 56.0582 




-35 45 32.5774 


3.357 


-2.7 


17 20 56.0580 




-35 45 32.5891 


6.453 


-3.0 


17 20 56.0581 




-35 45 32.5907 


11.608 


-3.4 


17 20 56.0579 




-35 45 32.5977 


16.950 


-3.7 


17 20 56.0579 




-35 45 32.5953 


13.072 


-4.0 


17 20 56.0579 




-35 45 32.5929 


10.177 


-4.4 


17 20 56.0579 




-35 45 32.5925 


13.265 


-4.7 


17 20 56.0577 




-35 45 32.5900 


16.513 


-5.0 


17 20 56.0578 




-35 45 32.5893 


20.479 


-5.4 


17 20 56.0578 




-35 45 32.5910 


28.570 


-5.7 


17 20 56.0579 




-35 45 32.5903 


38.006 


-6.0 


17 20 56.0578 




-35 45 32.5935 


38.297 


-6.3 


17 20 56.0578 




-35 45 32.5900 


28.243 


-6.7 


17 20 56.0577 




-35 45 32.5910 


19.718 


-7.0 


17 20 56.0578 




-35 45 32.5912 


22.008 


-7.3 


17 20 56.0578 




-35 45 32.5952 


31.096 


-7.7 


17 20 56.0577 




-35 45 32.5945 


35.008 


-8.0 


17 20 56.0579 




-35 45 32.5918 


28.781 


-8.3 


17 20 56.0579 




-35 45 32.5943 


18.181 


-8.6 


17 20 56.0580 




-35 45 32.5908 


10.865 


-9.0 


17 20 56.0582 




-35 45 32.5926 


7.395 


-9.3 


17 20 56.0582 




-35 45 32.5850 


4.048 


-9.6 


17 20 56.0582 




-35 45 32.5758 


1.460 


-10.0 


17 20 56.0600 




-35 45 32.6461 


0.521 


-15.6 


17 20 56.0558 




-35 45 32.4942 


0.536 


-15.9 


17 20 56.0566 




-35 45 32.5314 


0.608 


-16.2 


17 20 56.0558 




-35 45 32.5633 


0.515 


-21.5 


17 20 56.0600 




-35 45 32.3669 


0.375 


+5.2 


17 20 56.0585 




-35 45 32.4918 


0.432 


+4.9 


17 20 56.0597 




-35 45 32.5667 


0.562 


+4.5 


17 20 56.0587 




-35 45 32.5882 


0.681 


+4.2 


17 20 56.0579 




-35 45 32.5851 


1.028 


+3.9 


17 20 56.0582 




-35 45 32.5809 


1.426 


+3.5 


17 20 56.0586 




-35 45 32.5788 


1.805 


+3.2 


17 20 56.0579 




-35 45 32.6131 


1.706 


+2.9 


17 20 56.0580 




-35 45 32.5580 


1.577 


+2.5 


17 20 56.0581 




-35 45 32.5956 


1.228 


+2.2 


17 20 56.0584 




-35 45 32.5573 


0.963 


+1.9 


17 20 56.0583 




-35 45 32.5609 


0.652 


+1.6 


17 20 56.0613 




-35 45 32.4586 


0.509 


+0.9 


17 20 56.0556 




-35 45 32.5405 


0.566 


+0.6 


17 20 56.0580 




-35 45 32.5863 


0.547 


+0.2 


17 20 56.0595 




-35 45 32.5222 


0.460 


-0.1 


17 20 56.0603 




-35 45 32.4752 


0.459 


-0.4 


17 20 56.0583 




-35 45 32.5707 


0.450 


-0.7 


17 20 56.0600 




-35 45 32.5378 


0.579 


-1.1 


17 20 56.0565 




-35 45 32.6301 


0.636 


-1.4 


17 20 56.0588 




-35 45 32.5398 


0.573 


-1.7 


17 20 56.0581 




-35 45 32.6371 


0.745 


-2.1 


17 20 56.0587 




-35 45 32.5701 


1.514 


-2.4 


17 20 56.0582 




-35 45 32.5774 


3.370 


-2.7 


17 20 56.0580 




-35 45 32.5890 


6.416 


-3.0 


17 20 56.0581 




-35 45 32.5906 


11.555 


-3.4 


17 20 56.0579 




-35 45 32.5976 


16.873 


-3.7 


17 20 56.0579 




-35 45 32.5952 


13.029 
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TABLE 4 

44 GHz Class I Methanol masers detected in NGC 6334 1(N) 



Spot 


Peak position 


Peak intensity 


brightness 




Velocity (km s ') 




Number 


R.A. 


Dec 


(Jy beam"') 


(K) 


^peak 


^min 


Vmax 


1 


1 7 on aci 


-J J 4J LJ.Dl 


13.48 ±0.22 


16340 


-3.17 


-3.67 


-2.67 


2 


1 / ZU J J. 464 


'X'\ A'> 14 QO 

-J J 4j 14. yu 


5. 19 ±0.20 


6290 


-5.50 


-5.83 


-5.17 


3 




A'^ 18 98 
-JJ 4J io.Zo 


2.69 ± 0.20 


3260 


-4.17 


-4.34 


-3.84 


4 


i / ZU 3 J. 1 oy 


-JJ 4J iJ.4Z 


12.82±0.18 


15540 


-1.68 


-2.18 


-Hi. 14 


5 


1 / ZU Jt.WJZ 


AS IS Sf^ 

-JJ H-J iJ.JO 


0.82±0.17 


990 


-2.18 


-2.84 


-1.84 


6 


1 / ZU J^.UJJ 


^S AS 90 '^8 
-JJ H-J ZU. Jo 


3.01 ±0.18 


3650 


-4.50 


-5.17 


-4.17 


7 


1 7 7n 774 
i / ZU J4.Z Z^- 


^S AS 17 11 
-J J 4J iZ. i i 


1.10±0.16 


1330 


-7.16 


-7.66 


-6.49 


8 


1 7 7n 1^1 

i / ZU J4. J 


^S /I ^ 77 1 Q 

-J J 4j zz. iy 


10.71 ±0.17 


12980 


-6.83 


-7.82 


-6.16 


9 


1 7 7n Af\K 
i / ZU J^.^Uj 


"^S A^ 77 S7 
-J J 43 ZZ. J / 


4.39 ±0.17 


5320 


-8.49 


-9.98 


-7.99 


10 


1 7 70 Szl zl 1 Q 


^S AS 1 S SO 
-JJ 4J iJ.JU 


0.95 ±0.16 


1200 


-8.32 


-9.32 


-8.16 


11 


1 7 70 Szl zlS7 
i / ZU J^.^-J / 


^S AS 97 S7 
-J J 4J ZZ. J / 


83.28 ±0.17 


100900 


-7.16 


-8.16 


-6.16 


12 


M 0(\ '^A AlA 

1 1 ZU Jf.f / *+ 


-'^S AS 1 A Sf^ 


0.53±0.15 


640 


-1.51 


-1.84 


-0.85 


13 


1 7 90 '^A A7R 
1 / ZU .f / o 


"^S AS (\f\ f\A 


0.62 ±0.15 


750 


-1.35 


-1.84 


-1.01 


14 


1 7 on <A 

i / ZU J4. J JZ 


"^S AS 99 A(\ 
-Do 4J ZZ.40 


5.46 ±0.17 


6620 


-6.16 


-7.66 


-5.17 


15 


1 7 90 ^ifi9 

i / ZU Jf.UDZ 


'^S AS A'\ S9 

-JJ tj tJ.JZ 


5. 10 ±0.37 


6180 


-4.83 


-5.17 


-3.67 


16 


1 7 7n C/1 o-iQ 
i / ZU J^.ojo 


ac AS "^7 '^l 

-J J 4 J J / . J i 


6.63 ±0.26 


8040 


-4.83 


-5.33 


-4.67 


17 


1 7 70 79 
i / ZU J4. /ZJ 


^S AS 1 A 79 
-JJ 4J 14. ZZ 


18.83 ±0.15 


22820 


-5.33 


-6.00 


-3.67 


18 


1 7 70 7S1 
i / ZU J4. / o i 


■^S AS OQ fn 
-JJ 4J uy.D / 


0.81 ±0.15 


980 


-2.01 


-3.01 


-1.68 


19 


1 7 70 ^ /I Q7 1 
i / ZU j4.oZ i 


"IS /I ^ ] A Al 
-JJ 4j 14. 4j 


1.40 ±0.15 


1700 


-4.67 


-4.83 


-4.34 


20 


1 7 70 Szl 8A8 


"^S AS "^f^ 8A 
-JJ 4J JD.o4 


6.11 ±0.25 


7410 


-4.34 


-4.67 


-3.84 


21 


1 7 70 1 

i / ZU J4.00 1 


^S AS 11 ^^S 
-J J 4 J 1 1 .DJ 


0.74 ±0.14 


900 


-1.68 


-2.01 


-1.35 


22 


1 7 70 SA Q9 1 
i / ZU JH-.yz 1 


^S AS 1 '^ 88 
-J J H-J IJ.oo 


33.12±0.15 


40140 


-3.34 


-3.84 


-2.84 


23 


1 7 70 SA Q9S 

i / ZU JM-.^Zo 


■^s AS 1 '^ ^lO 

-JJ M-J IJ.OU 


12.19±0.15 


14770 


-5.17 


-6.49 


-4.83 


24 


1 7 70 SA Q'^^ 


AS 08 77 
-JJ 4j Uo.Z / 


3.89 ±0.14 


4710 


-4.67 


-5.00 


-4.34 


25 


1 7 on <A QftO 


o< AS 1 OA 
-JJ 4 J I J. 74 


78.63 ±0.15 


95300 


-4.34 


-4.67 


-4.00 


26 


1 7 90 ^A Q8A 
1 / ZU D4.yo4- 


"^S AS 1 OQ 
-JJ 4 J 1 j.Uo 


0.95 ±0.14 


1200 


-5.17 


-5.66 


-4.83 


27 


1 7 on <A QQl 
i / ZU j4.7yi 


"^S AS 9R R7 
-JJ 4 J Zo.o/ 


2.64 ±0.19 


3200 


-4.50 


-4.83 


-4.17 


28 


1 7 on c< nno 

i / ZU JJ.UU7 


'5< AS 1 A Oft 
-JJ 4j 14.UO 


35.40 ±0.15 


42900 


-2.67 


-3.67 


-1.84 


29 


1 7 70 OQO 


"^S AS 1 S SA 

-J J "4- J 1 J. JM- 


4.01 ±0.15 


4860 


-2.18 


-3.17 


-1.51 


30 


1 7 7f^ flQS 

1 / ZV/ JJ.V^^J 


zLS ^fi 78 
-JJ 4-J JO.Zo 


1.16±0.24 


1410 


-0.35 


-0.68 


-0.02 


31 


1 7 7f^ 1 OQ 

i / ZU J J. iuy 


^S AS 88 
-JJ 4J JD.oo 


1.86 ±0.25 


2250 


-7.16 


-7.49 


-6.66 


32 


1 7 70 1 

i / ZU J J. iOJ 


^S AS ^7 S^ 
-J J '-tJ J / .J J 


6.05 ± 0.26 


7330 


-3.67 


-4.00 


-2.84 


33 


1 7 70 1 Q'^ 
i / ZU J J. lyj 


^S AS 0'? 8 1 
-J J 4 J UJ.o 1 


1.02 ±0.1 5 


1240 


-0.85 


-2.01 


-0.02 


34 


1 7 70 708 
i / ZU JJ.ZUo 


^S AS ^7 SS 
-J J 4 J J / . Jo 


1 .34 ± 0.26 


1620 


-2.01 


-2.34 


-1.68 


35 


1 7 70 7^A 
i / ZU J J.ZO^ 


"^S AS 'if^ 8A 
-JJ 4J JD.o4 


4.79 ±0.25 


5810 


-3.84 


-4.34 


-2.67 


36 


1 7 70 7^1 
i / ZU J J.Zj 1 


-J J 4j JU. i o 


18.37±0.19 


22260 


-0.68 


-1.51 


-0.35 


37 


1 7 7n 'ifiA 
i / ZU JJ. juu 


S AS ^0 IS 
-J J 4J JU.lJ 


1.43 ±0.19 


1730 


-1.18 


-1.68 


-1.01 


38 


1 7 on cc aoQ 

1 / ZU JJ. JZo 


"^S AS 01 It^ 
-JJ 4 J Ul. /O 


1.42 ±0.15 


1720 


-1.84 


-2.18 


-1.18 


39 


1 7 on -ion 


"^S AS OS SS 
-JJ 4 J UJ.JJ 


0.60 ±0.15 


730 


-1.18 


-1.51 


-0.85 


40 


1 7 on cc 'if.-i 
i / ZU JJ. JD / 


'XC AS OS rtl 
-JJ 4 J UJ.Ol 


1.37±0.15 


1660 


-3.34 


-4.00 


-2.84 


41 


1 7 on cc -JQO 


'5C AS OS 1 
-JJ 4 J UJ.IO 


1.02 ±0.15 


1240 


-1.01 


-1.51 


-0.68 


42 


1 7 7n cc AAO 

i / ZU JJ.'-f4U 


•3C AS OQ 78 

-J J 4J uy. / o 


11.19±0.14 


13560 


-4.67 


-5.33 


-4.50 


43 


1 7 7n Af^f, 
i / ZU JJ.^OD 


^S HA A7 
-JJ 4j UD.DZ 


2.83 ±0.15 


3430 


-1.35 


-2.01 


-0.85 


44 


1 7 70 S'^ A71 
1 / ZU J J. 4 / 1 


■^S AS 00 iifx 
-JJ 4J UU.oO 


1.61 ±0.15 


1950 


-1-2.64 


-Hi. 64 


-H4.80 


45 


1 7 70 S'^ '^l 

i / ZU J J.J i J 


ac /1C (-\/:: 9c 
-JJ 4J UD.ZJ 


1.38±0.15 


1670 


-4.50 


-4.83 


-4.34 


46 


1 7 70 SI 

i / ZU J J.J i J 


"^S AS 07 AA 
-J J M-j u / .uM- 


1.90 ±0.15 


2300 


-5.00 


-5.50 


-4.67 


47 


1 7 70 S^S 

i / ZU J J. J Jo 


^S AS 00 ftO 

-J J "4 J UU.OZ 


1.63 ±0.15 


1980 


-4.17 


-4.34 


-3.84 


48 


1 7 70 '^A^ 
i / ZU J J. J4U 


^S AS OQ OS 

-JJ 4J uy.uj 


0.74 ±0.14 


900 


-5.33 


-5.50 


-5.17 


49 


1 7 7n c^o 

i / ZU J J. JJo 


^S AS 9f^ 79 
-J J 4J ZD. /Z 


3.32±0.18 


4020 


-3.34 


-3.67 


-3.01 


50 


17 90 "^Rl 

i / ZU JJ.Jol 


^S AS 07 S^ 
-JJ M-J UZ.JJ 


0.96 ±0.15 


1200 


-2.01 


-3.01 


-0.68 


51 


1 7 on <c n 

i / ZU JJ.DIU 


-JJ 4 J Z7.4J 


5.35±0.19 


6480 


-2.84 


-3.17 


-2.51 


52 


17 90 SS f\'\\ 


-'^S AS '^1 9"^ 

-JJ f J Jl.ZJ 


4.48 ±0.21 


5430 


-4.34 


-4.50 


-4.17 


53 


17 90 f\'^0 

1 / ZU JJ.OJZ 


"^S AS 01 fi,Si 

-JJ f J Ul.UO 


1.68±0.15 


2040 


-1.18 


-2.18 


-0.02 


54 


1 7 on f,iif. 

i / ZU JJ.UJU 


AS in rt9 
-JJ 4 J lU.OZ 


2.86±0.15 


3470 


-5.83 


-6.49 


-4.67 


55 


1 7 on c< ^:77 
i / ZU JJ.O / / 


'^s AS on 
-JJ 4 J uu.jy 


296.42 ±0.16 


359260 


-2.67 


-3.17 


-1.51 


56 


1 7 on c< f.Q'X 

i / ZU JJ.OoJ 


'^S AS 09 88 
-JJ 4 J UZ.oo 


1.02 ±0.15 


1240 


-7.16 


-7.66 


-6.66 


57 


1 7 7n cc 7nQ 
i / ZU J J. /Uo 


AS 1 1 17 

-J J 4 J 11.1/ 


1.42±0.15 


1720 


-3.51 


-3.67 


-3.17 


58 


1 7 70 771 
i / ZU J J. / / 1 


ac A'Z 77 7S 
-JJ 4 J Z / . / J 


268.96 ±0.19 


325980 


-4.50 


-5.17 


-3.84 


59 


1 7 70 8 1 

i / ZU J J.o iU 


■^S AS "^1 07 
-J J 4J J i .U / 


2.12±0.21 


2570 


-5.00 


-5.66 


-4.83 


60 


1 7 7n 

i / ZU JJ.ojy 


"5 S /I ^ ^1 7/1 
-J J 43 J i . / 4 


3.45 ±0.21 


4180 


-5.33 


-6.00 


-4.83 


61 


1 7 70 8SQ 

i / ZU J j.ojy 


-ic /ic 90 97 
-J J 4J Zo.Z / 


7.12±0.19 


8629 


-3.84 


-4.34 


-3.51 


62 


1 7 70 noA 

i / ZU JD.UUt- 


"^S AS ^'^ OA 

-J J M-J JJ.UM- 


4.53 ±0.23 


5490 


-3.84 


-4.50 


-3.17 


63 


1 7 70 Sf^ OOA 

i / ZU JU.IAJH- 


"^S AS 9'^ A9 

-JJ M-J ZJ.M-Z 


1.98 ±0.1 8 


2400 


-3.34 


-3.67 


-3.17 


64 


1 7 70 flQS 
i / ZU JDAJVJ 


99 /II 
-JJ 4J ZZ. 41 


4.12±0.17 


4990 


-4.67 


-5.00 


-4.00 


65 


1 7 70 1 98 
i / ZU JO. IZo 


/IC 90 CA 
-JJ 4J ZJ.JU 


28.48±0.18 


34520 


-5.00 


-5.66 


-4.50 


66 


1 7 90 '\(^ 1 
i / ZU DO. 104 


OC A< 00 80 
-JJ 4 J ZZ. oU 


5.23 ±0.18 


6340 


-7.66 


-8.65 


-7.32 


67 


1 7 on 1 Q< 
1 / ZU JO. loj 


-ic AS O'X 8S 
-JJ 4 J ZJ.oJ 


10.17±0.18 


12330 


-7.82 


-8.82 


-6.49 


68 


17 20 56.196 


-35 45 32.22 


9.45 ±0.23 


11450 


-4.17 


-4.50 


-3.67 


69 


17 20 55.264 


-35 45 36.84 


4.80 ±0.25 


5820 


-5.33 


-5.66 


-5.00 


70 


17 20 56.285 


-35 45 19.20 


2.14±0.17 


2590 


-3.51 


-4.34 


-3.17 


71 


17 20 56.296 


-35 45 29.73 


43.26 ±0.22 


52430 


-5.00 


-5.83 


-4.50 


72 


17 20 56.342 


-35 45 20.41 


1.74±0.18 


2110 


-3.51 


-4.17 


-3.17 


73 


17 20 56.307 


-35 45 26.75 


0.90 ±0.20 


1100 


-2.18 


-2.51 


-2.01 


74 


17 20 56.341 


-35 45 30.68 


28.14 ±0.23 


34110 


-4.00 


-4.50 


-3.17 


75 


17 20 56.414 


-35 45 31.72 


171.47±0.24 


207820 


-3.34 


-3.67 


-3.01 


76 


17 20 56.423 


-35 45 29.69 


9.90 ±0.23 


12000 


-4.00 


-6.00 


-3.67 


77 


17 20 56.488 


-35 45 30.49 


5.38 ±0.24 


6520 


-2.51 


-3.67 


-2.01 


78 


17 20 56.708 


-35 45 32.61 


13.34 ±0.28 


16170 


-2.51 


-3.01 


-1.18 


79 


17 20 57.028 


-35 45 35.92 


2.02 ±0.37 


2450 


-5.33 


-5.66 


-4.34 


80 


17 20 56.903 


-35 45 24.93 


2.06 ±0.24 


2500 


-3.51 


-3.84 


-3.17 
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TABLE 5 

1.3 MILLIMETER SPECTRAL LINES DETECTED TOWARD SMAl 



Species 


Transition 


Frequency 


(cm 

'-'lower v^A'A / 


Brightness 


H9S 


22,0 -2i,i 


216.710435 


51.1 


13 


HCOOCH3-E 


184,15-174,14 


216.830112 


66.2 


6.2 


HLUULH3-A 


182,16-172,15 


216.838646 


66.2 


4.8 


C2H3CN 


^ IT Z,Zj t t 

^^Z,ZZ ^-^,z,zl 


216.936717 


85.9 


0.8 


CH3OH-E 


^+1,4 4+2,2 


216.945600 


31.6 


11 


HCOOCH3-E 




216.964789 


70.3 


10 




20l,20 — 1^1,19 


21 fi 965991 


70.3 


Il3 


HCOOCH3-E 


202,19 -19o,i9 


216.966247 


70.3 


J 


HCOOCH3-A 


20o,20-l%,19 


216.967334 


70.3 


11 


Unidentified 


217.020 




2.3 


SiO 


5-4 


217.104980 


14.5 


2.8 


CH3OCH3 


36o,36,l -36l,35,l 


217.168120 


436.6 


1.2 


CH3OCH3 


224,19,2-223,20,2 


217.189667 


168.9 


1 


CH3OCH3 




217.189667 


168.9 


!> 2.9 


CH3OCH3 


2.Z,A 101 — 2^2.1 on 1 

4, ly, 1 -^^jiZU, 1 


217.191424 


168.9 


4.8 


CH3OCH3 


Z 2 A inn — 2.2.1 inn 


217.193170 


168.9 


3.6 


DCN 


3-2'^ 


217.238630 


7.2 


5.7 


CH3OH-A- (v, = 1) 


61,5 -72.6 


217.299202 


252.6 


8.4 


Unidentified 


217.312 




2.2 


''CH3OH-A+ 


102,8-93,7 


217.399550 


105.6 


2.4 


Unidentified 


217.418 




1.9 


CH3OH-A+ (v, = 1) 


1 ii; m — 1 lie 1 1 


217.642859 


511.0 


3.3 


CH3OH-E 


20+1^19 ~ 20+0^20 


217 886390 


346.1 


5.0 


Unidentified 


217.959 




1.9 


oc"s 


18—17 


718 198998 


61 9 


3 7 




3o,3 ~2o.2 


718 979197 


7 3 


8 4 


1. 1 v_y V_>' V 1 ±^ I -> 


A'1,16 ^^^1,15 


218 2808^5 


62.0 


4.4 




173,14 — 163,13 




fO ft 






Z4 — 23 




o3.o 


21 


C2H5CN 


243,21-233,20 


zlo.joyy /U 


yu.u 


1.1 




'++2,2- J+1,2 


1 c /i/inH/in 
Z16.44UU4U 


24.3 


23 




^^5,19 -^^5,18 


Zlo.4jlzyo 


11/./ 


2. / 




235.18-225,17 




117 7 


2. J 




iui,9-yi,8 




K A 
3j.U 


^ /I 
3.4 


H2CO 


32,2 22,1 


Z i 0.'+ / JXJjZ, 


40 4 


7 7 


CH3OCH3 


233,21,1-232,22,1 


218.492412 


176.1 


4.7 


C2H3CN 


234,20-224,19 


218.573646 


104.2 


1.7 


C2H3CN 


233,21 -223,20 


218.585072 


93.7 


1.6 


CN 


2o,2,3 - lo,l,2 


226.659558 


3.8 


-3.4 


CN 


2o,2,i - lo,i,i 


226.663693 


3.8 


-0.9 


CN 


2o,2,2 - lo.l.l 


226.679311 


3.8 


-4.9 


HCOOCH3-E 


202,18-192,18 


226.713060 


76.0 


5.8 


HCOOCH3-A 


202,19-192,18 


226.718688 


76.0 


5.3 


HCOOCH3-E 


203,18-193,17 


226.773130 


76.0 


5.2 


HCOOCH3-A 


20l,l9-19i,i8 


226.778707 


76.0 


6.0 


CN 


2o,3,3 - lo,2,2 


226.874191 


3.8 


1 




2o,3,4- lo,2,3 


ZZO.o / 4 / cS i 


■? s 
3.0 


\ -6.9 
J 


CN 


2o,3,2 — lo,2,l 


Z.Z.U . / J y\j 


3.8 


-6.1 


CN 


2o,3,2 — lo,2,2 


226.887420 


3.8 


-2.2 


CN 


2o,3,3 — lo,2,3 


226.892128 


3.8 


-1.5 


HCOOCH3-E 


194,16-184,15 


227.019490 


73.5 


4.8 


HCOOCH3-A 


192,17-182,16 


227.028070 


73.4 


5.0 


CH3OH-E 


21+1,20-21+0,21 


227.094601 


379.6 


6.9 


Unidentified 


227.126 




2.2 


Unidentified 




22 /.I /O 




2.0 


Unidentified 




22 / .221 




1 .& 


CH3OH-E 


1 2 1 n — 11 

^■^^i,iz -■- -"-H-Zjy 


227.229600 


122.0 


2.2 


g-CzHjOH 


133,10.1 - 123,9.1 


227.294752 


95.7 


1.8 


13c34s 


5-4 


227.300506 


15.16 


3.9 


HC3N 


25 - 24'' 


227.418905 


91.03 


17 


HCOOCH3-E 


21l,21 -20l,20 


227.560955 


77.5 


12 


1. 1 V_^V_^V_^ 1 ±3 L ' 


212,20 — 202,19 


227 56175^ 


77.5 


I12 


HCOOCH3-A 


21l,21 -20i,20 


227.561944 


77.5 


/ 


HCOOCH3-A 


21o,21 -20o,20 


227.562740 


77.5 


10 


t-C2H50H 


185,13,2-184,14,2 


227.606079 


114.3 


l.J 


NH2CHO 


1 lo 11 ~ IOq 10 


227.606176 


38.5 




LH3UCH3 


265,21,0-264,22,0 


227.047849 


239.7 




CH3OCH3 


265,21,1 -264,22,1 


227.647984 


239.7 


}« 


CH3OCH3 


265,21,2-264,22,2 


227.648120 


239.7 




CH3OCH3 


243,22,2 - 242,23,2 


227.654386 


190.9 




CH3OCH3 


243,22.3 -242.23,3 


227.654386 


190.9 


} 4.6 


CH3OCH3 


243,22,1 -242,23,1 


227.657025 


190.9 


4.9 


CH3OCH3 


243,22,0-242,23,0 


227.659664 


190.9 


2.0 


C2H5CN 


253,22-243,21 


227.780972 


97.3 


6.1 


CH3OH-A+ 


I61, 16- 152,13 


227.814650 


219.8 


9.5 


g-CzHjOH 


13l,12,l - 12l,ll,l 


227.891911 


89.7 


0.9 




■74 O'X, .„ 


■707 QflfifiS^ 


141 S 
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TABLE 6 

Derived 1.3 millimeter Continuum Properties 



Source 




^gas 


Nh2 X 10^3 


nsj X 10'' 




(K) 


(Mo) 


(cm-2) 


(cm-3) 


SMAl 


155 - 165 


16- 14 


7.7 - 6.8 


1.3 - 1.2 


SMA2 


135 - 155 


7-6 


6.4 - 5.7 


1.6- 1.4 


SMA3 


20-50 


74-23 


40-12 


7-2 


SMA4 


20-50 


52- 16 


30-9 


6-2 


SMA5 


20-50 


28-9 


18-6 


4-1 


SMA6 


25-50 


32-13 


58-17 


16-5 


SMA7 


20-50 


22-7 


8-3 


1-0.4 



